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Abstract.  The  enhancement  in  diffusive  transport  of  passive  tracer  particles  by  incompressible,  turbulent  flow  fields  is 
a  challenging  problem  with  theoretical  and  practical  importance  in  many  areas  of  science  and  engineering,  ranging  from  the 
transport  of  mass,  heat,  and  pollutants  in  geophysical  flows  to  sea  ice  dynamics  and  turbulent  combustion.  The  long  time, 
large  scale  behavior  of  such  systems  is  equivalent  to  an  enhanced  diffusive  process  with  an  effective  diffusivity  tensor  D*.  Two 
different  formulations  of  integral  representations  for  D*  were  developed  for  the  case  of  time-independent  fluid  velocity  fields, 
involving  spectral  measures  of  bounded  self-adjoint  operators  acting  on  vector  fields  and  scalar  fields,  respectively.  Here,  we 
extend  both  of  these  approaches  to  the  case  of  space-time  periodic  velocity  fields,  with  possibly  chaotic  dynamics,  providing 
rigorous  integral  representations  for  D*  involving  spectral  measures  of  unbounded  self-adjoint  operators.  We  prove  that  the 
different  formulations  are  equivalent.  Their  correspondence  follows  from  a  one-to-one  isometry  between  the  underlying  Hilbert 
spaces.  We  also  develop  a  Fourier  method  for  computing  D*,  which  captures  the  phenomenon  of  residual  diffusion  related  to 
Lagrangian  chaos  of  a  model  flow.  This  is  reflected  in  the  spectral  measure  by  a  concentration  of  mass  near  the  spectral  origin. 
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1.  Introduction.  The  long  time,  large  scale  motion  of  diffusing  particles  or  tracers  being  advected  by 
an  incompressible  flow  field  is  equivalent  to  an  enhanced  diffusive  process  [94]  with  an  effective  diffusivity 
tensor  D*.  Describing  the  associated  transport  properties  is  a  challenging  problem  with  a  broad  range  of 
scientific  and  engineering  applications,  such  as  stellar  convection  [47,  82,  22,  23,  21],  turbulent  combustion  [3, 
17,  93,  99,  78,  101],  and  solute  transport  in  porous  media  [14,  15,  98,  41,  48,  51,  49].  Time-dependent  flows 
can  have  fluid  velocity  fields  with  chaotic  dynamics,  which  gives  rise  to  turbulence  that  greatly  enhances  the 
mixing,  dispersion,  and  large  scale  transport  of  diffusing  scalars. 

In  the  climate  system  [27,  40],  turbulence  plays  a  key  role  in  transporting  mass,  heat,  momentum, 
energy,  and  salt  in  geophysical  flows  [67].  Turbulence  enhances  the  dispersion  of  atmospheric  gases  [29]  such 
as  ozone  [43,  79,  80,  81]  and  pollutants  [26,  11,  86],  as  well  as  atmosphere-ocean  transfers  of  carbon  dioxide 
and  other  climatically  important  trace  gas  fluxes  [103,  8] .  Longitudinal  dispersion  of  passive  scalars  in  oceanic 
flows  can  be  enhanced  by  horizontal  turbulence  due  to  shearing  of  tidal  currents,  wind  drift,  or  waves  [102, 
50,  19].  Chaotic  motion  of  time-dependent  fluid  velocity  fields  cause  instabilities  in  large  scale  ocean  currents, 
generating  geostrophic  eddies  [33]  which  dominate  the  kinetic  energy  of  the  ocean  [34],  Geostrophic  eddies 
greatly  enhance  [33]  the  meridional  mixing  of  heat,  carbon  and  other  climatically  important  tracers,  typically 
more  than  one  order  of  magnitude  greater  than  the  mean  flow  of  the  ocean  [89].  Eddies  also  impact  heat 
and  salt  budgets  through  lateral  fluxes  and  can  extend  the  area  of  high  biological  productivity  offshore  by 
both  eddy  chlorophyll  advection  and  eddy  nutrient  pumping  [24]. 

In  sea  ice  dynamics,  where  the  ice  cover  couples  the  atmosphere  to  the  polar  oceans  [96],  the  transport 
of  sea  ice  can  also  be  enhanced  by  eddie  fluxes  and  large  scale  coherent  structures  in  the  ocean  [97,  54].  In 
sea  ice  thermodynamics,  the  temperature  field  of  the  atmosphere  is  coupled  to  that  of  the  ocean  through  sea 
ice,  which  is  a  composite  of  pure  ice  with  brine  inclusions  whose  volume  fraction  and  connectedness  depend 
strongly  on  temperature.  Convective  brine  flow  through  the  porous  microstructure  can  enhance  thermal 
transport  through  the  sea  ice  layer  [55,  100,  52]. 

It  has  been  noted  in  various  geophysical  contexts  [80,  81]  that  eddy-induced,  skew-diffusive  tracer  fluxes, 
directed  normal  to  the  tracer  gradient  [64],  are  generally  equivalent  to  antisymmetric  components  in  the 
effective  diffusivity  tensor  D*,  while  the  symmetric  part  of  D*  represents  irreversible  diffusive  effects  [83, 
87,  39]  directed  down  the  tracer  gradient.  The  mixing  of  eddy  fluxes  is  typically  non-divergent  and  unable 
to  affect  the  evolution  of  the  mean  flow  [64],  and  do  not  alter  the  tracer  moments  [39].  In  this  sense,  the 
mixing  is  non-dissipative,  reversible,  and  sometimes  referred  to  as  stirring  [28,  39].  Both  numerical  and 
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observational  studies  of  scalar  transport  have  suggested  that  tracers  are  advected  over  large  scales  by  a  fluid 
velocity  field  that  is  different  from  the  mean  flow  [75].  This  suggests  that  the  effective  diffusivity  tensor  D* 
should  be  spatially  and  possibly  also  temporally  inhomogeneous  [75]. 

Due  to  the  computational  intensity  of  detailed  climate  models  [40,  96,  70],  a  coarse  resolution  is  necessary 
in  numerical  simulations  and  parameterization  is  used  to  help  resolve  sub-grid  processes,  such  as  turbulent 
entrainment-mixing  processes  in  clouds  [53] ,  atmospheric  boundary  layer  turbulence  [20] ,  atmosphere-surface 
exchange  over  the  sea  [30]  and  sea  ice  [88,  1,  2,  95],  and  eddies  in  the  ocean  [60,  37].  In  this  way,  only  the 
effective  or  averaged  behavior  of  these  sub-grid  processes  are  included  in  the  models.  Here,  we  study  the 
effective  behavior  of  advection  enhanced  diffusion  by  time-dependent  fluid  velocity  fields,  with  possibly 
chaotic  dynamics,  which  gives  rise  to  such  a  parameterization,  namely,  the  effective  diffusivity  tensor  D*  of 
the  flow. 

In  recent  decades,  a  broad  range  of  mathematical  techniques  have  been  developed  which  reduce  the 
analysis  of  enhanced  diffusive  transport  by  complex  fluid  velocity  fields  with  rapidly  varying  structures  in 
both  space  and  time,  to  solving  averaged  or  homogenized  equations  that  do  not  have  rapidly  varying  data,  and 
involve  an  effective  parameter  [72,  61,  9,  16,  31,  71,  32,  59,  75,  77,  25,  42,  44,  57,  58,  101].  Motivated  by  [73], 
it  was  shown  [61]  that  the  homogenized  behavior  of  the  advection-diffusion  equation  with  a  random,  time- 
independent,  incompressible,  mean-zero  fluid  velocity  field,  is  given  by  an  inhomogeneous  diffusion  equation 
involving  the  symmetric  part  of  an  effective  diffusivity  tensor  D*.  Moreover,  a  rigorous  representation  of 
D*  was  given  in  terms  of  an  auxiliary  “cell  problem”  involving  a  curl- free  random  field  [61].  We  stress  that 
the  effective  diffusivity  tensor  D*  is  not  symmetric  in  general.  However,  only  its  symmetric  part  appears  in 
the  homogenized  equation  for  this  formulation  of  the  effective  transport  properties  of  advection  enhanced 
diffusion  [61]. 

The  incompressibility  condition  of  the  time-independent  fluid  velocity  field  was  used  [4,  5]  to  transform 
the  cell  problem  in  [61]  into  the  quasi-static  limit  of  Maxwell’s  equations  [46,  38],  which  describe  the  transport 
properties  of  an  electromagnetic  wave  in  a  composite  material  [66].  The  analytic  continuation  method  for 
representing  transport  in  composites  [38]  provides  Stieltjes  integral  representations  for  the  bulk  transport 
coefficients  of  composite  media,  such  as  electrical  conductivity  and  permittivity,  magnetic  permeability,  and 
thermal  conductivity  [66].  This  method  is  based  on  the  spectral  theorem  [92,  84]  and  a  resolvent  formula  for, 
say,  the  electric  field,  involving  a  random  self-adjoint  operator  [38,  69]  or  matrix  [68].  Based  on  [38],  the  cell 
problem  was  transformed  into  a  resolvent  formula  involving  a  bounded  self-adjoint  operator,  acting  on  the 
Hilbert  space  of  curl- free  random  vector  fields  [4,  5].  This,  in  turn,  led  to  a  Stieltjes  integral  representation 
for  the  symmetric  part  of  the  effective  diffusivity  tensor  D*,  involving  the  Peclet  number  Pe  of  the  flow  and 
a  spectral  measure  of  the  operator  [4,  5] .  A  key  feature  of  the  method  is  that  parameter  information  in  Pe 
is  separated  from  the  complicated  geometry  of  the  time- independent  flow,  which  is  encoded  in  the  measure 
fi.  This  property  led  to  rigorous  bounds  [5]  for  the  diagonal  components  of  D*.  Bounds  for  D*  can  also  be 
obtained  using  variational  methods  [5,  31,  71,  32]. 

The  mathematical  framework  developed  in  [61]  was  adapted  [75,  62,  57]  to  the  case  of  a  periodic, 
time-dependent,  incompressible  fluid  velocity  field  with  non-zero  mean.  The  velocity  field  was  modeled 
as  a  superposition  of  a  large-scale  mean  flow  with  small-scale  periodically  oscillating  fluctuations.  It  was 
shown  [75]  that,  depending  on  the  strength  of  the  fluctuations  relative  to  the  mean  flow,  the  effective 
diffusivity  tensor  D*  can  be  constant  or  a  function  of  both  space  and  time.  When  D*  is  constant,  only  its 
symmetric  part  appears  in  the  homogenized  equation  as  an  enhancement  in  the  diffusivity.  However,  when  D* 
is  a  function  of  space  and  time,  its  antisymmetric  part  also  plays  a  key  role  in  the  homogenized  equation.  In 
particular,  the  symmetric  part  of  D*  appears  as  an  enhancement  in  the  diffusivity,  while  both  the  symmetric 
and  antisymmetric  parts  of  D*  contribute  to  an  effective  drift  in  the  homogenized  equation.  The  effective 
drift  due  to  the  antisymmetric  part  is  purely  sinusoidal,  thus  divergence- free  [75].  This  is  consistent  with 
what  has  been  observed  in  geophysical  flows  in  the  climate  system,  as  discussed  above.  In  [62],  this  result 
was  extended  to  weakly  compressible,  anelastic  fluid  velocity  fields. 

Based  on  [14],  the  cell  problem  discussed  in  [75]  was  transformed  into  a  resolvent  formula  involving  a 
self-adjoint  operator,  acting  on  the  Sobolev  space  [63,  35]  of  spatially  periodic  scalar  fields,  which  is  also 
a  Hilbert  space.  In  the  case  where  the  mean  flow  and  periodic  fluctuations  are  time- independent,  the  self- 
adjoint  operator  is  compact  [14],  hence  bounded  [90].  This  led  to  a  discrete  Stieltjes  integral  representation 
for  the  antisymmetric  part  of  D*,  involving  the  Peclet  number  of  the  steady  flow  and  a  spectral  measure  of 
the  operator. 
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Here,  we  generalize  both  of  the  approaches  described  in  [4,  5]  and  [75]  to  the  case  of  a  periodic,  time- 
dependent  fluid  velocity  field,  allowing  for  chaotic  dynamics.  In  particular,  for  each  approach,  we  provide 
Stieltjes  integral  representations  for  both  the  symmetric  and  antisymmetric  parts  of  the  effective  diffusiv- 
ity  tensor  D*,  involving  a  spectral  measure  of  a  self-adjoint  operator.  In  this  time-dependent  setting,  the 
underlying  operator  becomes  unbounded.  The  spectral  theory  of  unbounded  operators  is  more  subtle  and 
technically  challenging  than  that  of  bounded  operators.  For  example,  the  domain  of  an  unbounded  oper¬ 
ator  and  its  adjoint  plays  a  central  role  in  the  spectral  characterization  of  the  operator.  Neglecting  such 
important  mathematical  details,  the  Stieltjes  integral  representation  for  D*  given  in  [4,  5]  was  extended 
to  the  time-dependent  setting  in  [6].  Here,  we  provide  a  mathematically  rigorous  formulation  of  Stieltjes 
integral  representations  for  D*  in  the  time-dependent,  unbounded  operator  setting.  We  prove  that  the  two 
approaches  described  in  [4,  5]  and  [75]  are  equivalent  in  this  setting,  and  that  their  correspondence  follows 
from  a  one-to-one  isometry  between  the  underlying  Hilbert  spaces.  We  also  establish  a  direct  correspondence 
between  the  effective  parameter  problem  for  D*  and  that  arising  in  the  analytic  continuation  method  for 
composite  media. 

Analytical  calculations  of  the  spectral  measure  underlying  the  effective  diffusivity  tensor  D*  have  been 
obtained  only  for  a  handful  of  simple  models  of  periodic  fluid  velocity  fields  to  date  such  as  shear  flows.  We 
develop  a  Fourier  method  for  the  computation  of  D*.  In  particular,  we  compute  the  effective  properties  for 
the  following  space-time  periodic  flow  in  two  spatial  dimensions,  with  x  =  (x,y), 

(1.1)  u(t,x)  =  (cosy, cosa;)  +  Ocost  (sin y, sin x),  0e(0, 1]. 

The  steady  part  (cosy,  cos x)  of  the  flow  is  subject  to  a  time-periodic  perturbation  that  causes  transition  to 
Lagrangian  chaos  [16,  104].  In  the  advection  dominated  regime,  we  shall  compare  our  computations  of  the 
effective  diffusivity  for  the  steady  6  =  0  and  dynamic  9=1  settings. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2  the  theory  of  homogenization  for  the  advection- 
diffusion  equation  for  space-time  periodic  flows  is  reviewed.  Novel  Stieltjes  integral  representations  for  the 
effective  diffusivity  tensor  are  also  obtained  for  a  large  class  of  fluid  velocity  fields,  involving  a  spectral 
measure  of  an  unbounded  self-adjoint  operator.  In  Section  3  we  use  Fourier  methods  to  transform  the 
eigenvalue  problem  for  the  operator  into  an  infinite  set  of  algebraic  equations.  This  provides  a  rigorous 
mathematical  foundation  for  the  computation  of  effective  diffusivities  for  such  flows.  This  framework  is 
employed  in  Section  4  to  compute  the  discrete  component  of  the  spectral  measure  associated  with  the  fluid 
velocity  field  in  (1.1)  for  time-independent  flow  9  =  0  and  time-dependent  flow  9  =  1.  Our  computations 
highlight  that  the  behavior  of  the  measure  near  the  spectral  origin  governs  the  behavior  of  the  effective 
diffusivity  in  the  advection  dominated  regime  of  small  molecular  diffusion.  In  particular,  we  demonstrate 
that  for  9  =  0  there  is  a  spectral  gap  in  the  measure  near  a  limit  point  at  the  spectral  origin,  giving  rise  to 
the  known  vanishing  asymptotic  behavior  of  2D  cell  flows  [31,  71].  However  in  the  time  dependent  setting, 
a  strong  concentration  of  measure  mass  near  the  spectral  origin  gives  rise  to  the  phenomenon  of  residual 
diffusivity  in  the  limit  of  vanishing  molecular  diffusion. 

The  proofs  of  the  integral  representation  formulas  for  the  effective  diffusivity,  over  both  vector  and 
scalar  fields,  as  well  as  background  material,  are  in  the  appendices.  The  spectral  theory  of  unbounded  self- 
adjoint  operators  in  Hilbert  space  is  reviewed  in  Sections  A  and  B.  Two  mathematical  formulations  of  the 
effective  parameter  problem  for  advection  enhanced  diffusion  are  presented  in  Section  C,  leading  to  novel 
integral  representations  for  the  symmetric  and  antisymmetric  components  of  the  effective  diffusivity  tensor. 
In  Section  D  we  use  powerful  methods  of  functional  analysis  to  prove  that  the  two  approaches  are  equivalent, 
which  follows  from  a  one-to-one  isometry  between  the  associated  Hilbert  spaces.  In  Section  D.l  we  derive 
an  explicit  formula  for  the  discrete  component  of  the  spectral  measure,  which  is  employed  in  our  numerical 
computations. 

2.  Effective  transport  by  advective-diffusion.  The  density  6  of  a  cloud  of  passive  tracer  particles 
diffusing  along  with  molecular  diffusivity  e  and  being  adverted  by  an  incompressible  velocity  field  u  satisfies 
the  advection-diffusion  equation 

(2.1)  dtc/)(t,  x)  =  u(t,  x)  •  V0(f,  x)  +  eA(j)(t,  x),  ^>(0,  x)  =  4>o(x), 

for  t  >  0  and  x  £  Rd.  Here,  the  initial  density  (f>o(x)  and  the  fluid  velocity  field  u  are  assumed  given,  and 
u  satisfies  V  •  u  =  0.  In  equation  (2.1),  the  molecular  diffusion  constant  e  >  0,  d  is  the  spatial  dimension  of 
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the  system,  dt  denotes  partial  differentiation  with  respect  to  time  t,  and  A  =  V  •  V  =  V2  is  the  Laplacian. 
Moreover,  i/>  •  ip  =  xpTTp,  xj}T  denotes  transposition  of  the  vector  xf>,  and  Tp  denotes  component- wise  complex 
conjugation,  with  x/j  •  xf)  =  |i/>|2.  We  stress  that  all  quantities  considered  in  this  section  are  real-valued. 

We  consider  enhanced  diffusive  transport  by  a  periodic  fluid  velocity  field  and  non-dimensionalize  equa¬ 
tion  (2.1)  as  follows.  Let  d  and  T  be  typical  length  and  time  scales  associated  with  the  problem  of  in¬ 
terest.  Mapping  to  the  non-dimensional  variables  t  i— >■  t/T  and  x  K >-  x/d:  one  finds  that  <f>  satisfies  the 
advection-diffusion  equation  in  (2.1)  with  a  non-dimensional  molecular  diffusivity  e  i— >  Te/d2  and  veloc¬ 
ity  field  u  i — )■  T  u /  d.  There  are  several  different  non-dimensionalizations  possible  for  the  advection-diffusion 
equation.  A  detailed  discussion  of  various  non-dimensionalizations  involving  the  Strouhal  number,  the  Peclet 
number,  and  the  periodic  Peclet  number  is  given  in  [62,  57].  Here,  we  focus  on  the  long  time,  large  scale 
transport  characteristics  of  equation  (2.1)  as  a  function  of  e.  To  this  end,  we  simply  take  T  to  be  the 
temporal  periodicity  of  the  velocity  field  u  and  assume  that  the  spatial  periodicity  of  u  is  d  in  all  spatial 
dimensions,  i.e. , 

(2.2)  u(t +  T,x)  =  u(t,x),  u(t,x  +  dej)  =  xi(t,x),  j  =  l,...,d, 

where  e7  is  a  standard  basis  vector  in  the  jth  direction. 

The  long  time,  large  scale  dispersion  of  diffusing  tracer  particles  being  advected  by  an  incompressible 
fluid  velocity  field  is  equivalent  to  an  enhanced  diffusive  process  [94]  with  an  effective  diffusivity  tensor  D*. 
In  recent  decades,  methods  of  homogenization  theory  [61,  31,  71,  57]  have  been  used  to  provide  an  explicit 
representation  for  D*.  In  particular,  these  methods  have  demonstrated  that  the  averaged  or  homogenized 
behavior  of  the  advection-diffusion  equation  in  (2.1),  with  space-time  periodic  velocity  field  «,  is  determined 
by  a  diffusion  equation  involving  an  averaged  scalar  density  </>  and  an  effective  diffusivity  tensor  D*  [57] 

(2.3)  3^(t,®)  =  V.[D*V^(t,*)],  0(0,  cc)  =  4>o(x). 

Equation  (2.3)  follows  from  the  assumption  that  the  initial  tracer  density  (f> o  varies  slowly  relative  to  the 
variations  of  the  fluid  velocity  field  u  [61,  32,  57].  This  information  is  incorporated  into  equation  (2.1)  by 
introducing  a  small  dimensionless  parameter  S  <C  1  and  writing  [61,  32,  57] 

(2.4)  c/)(0,x)  =  0o(Sx). 

Anticipating  that  will  have  diffusive  dynamics  as  t  — >  oo,  space  and  time  are  rescaled  according  to  the 
standard  diffusive  relation 

(2.5)  £  =  x/5,  t  =  t/51 ,  7  =  2. 

The  rescaled  form  of  equation  (2.1)  is  given  by  [57] 

(2.6)  dt4>s(t,  x)  =  5~1u(t/52,  x/5)  •  Vc/)6(t,  x)  +  sAcj)S(t,  x),  </>(0,  x)  =  4>o(x), 

where  we  have  denoted  x)  =  <j>(t/5 2,  x/5).  The  convergence  of  4>s  to  0  can  be  rigorously  established  in 
the  following  sense  [57] 

(2.7)  lim  sup  sup  |</>5(t,  x)  —  x)\  =  0, 

0<t<t0  a;GRd 

for  every  finite  to  >  0,  provided  that  <j>o  and  u  obey  some  mild  smoothness  and  boundedness  conditions,  and 
that  u  is  mean-zero. 

For  fixed  0  <  5  <C  1,  an  explicit  representation  of  the  effective  diffusivity  tensor  D*  is  given  in  terms  of 
the  (unique)  mean  zero,  space-time  periodic  solution  \j  of  the  following  cell  problem  [16,  57], 

(2-8)  dTXj(T,0  -  eAzXj(T,£)  -  u(r,£)  ■  V^Xj(t,0  =  Uj(r,£), 

where  the  subscript  £  in  and  Vj  indicates  that  differentiation  is  with  respect  to  the  fast  variable  £ 
defined  in  equation  (2.5).  Specifically,  the  components  D*fc,  j,k  =  1  of  the  matrix  D*  are  given 

by  [61,  31,  71,  57] 


(2.9) 


Djk  =  eSjk  +  (UjXk), 


Spectral  analysis  of  advection-diffusion  by  dynamic  flows 


5 


where  5jk  is  the  Kronecker  delta  and  u3  is  the  jth  component  of  the  vector  u.  The  averaging  (•)  in  (2.9) 
is  with  respect  to  the  fast  variables  defined  in  equation  (2.5).  More  specifically,  consider  the  bounded  sets 
T  Cl  and  V  C  Rd,  with  r  £  T  and  £  £  V,  which  define  the  space-time  period  cell  ((d+  l)-torus)  T  x  V.  In 
the  case  of  a  time-dependent  fluid  velocity  field,  (•)  denotes  space-time  averaging  over  Tx  V.  In  the  special 
case  of  a  time- independent  fluid  velocity  field,  the  function  Xj  is  time-independent  and  satisfies  equation  (2.8) 
with  dTXj  =  0,  and  (•)  in  (2.9)  denotes  spatial  averaging  over  V  [31,  71,  57].  Since  the  remainder  of  the 
analysis  involves  only  the  fast  variables,  for  notational  simplicity,  we  will  drop  the  subscripts  £  displayed  in 
equation  (2.8). 

In  general,  the  effective  diffusivity  tensor  D*  has  a  symmetric  S*  and  antisymmetric  A*  part  defined  by 

(2.10)  D*  =  S*+A*,  S*  =  i(D*  +  [D*]T),  A*  =  i  (D*  —  [D*] T) , 

where  [D*] 1  denotes  transposition  of  the  matrix  D*.  Denote  by  S*k  and  A*k,  j,  k  =  1, . . . ,  d,  the  components 
of  S*  and  A*  in  (2.10).  In  Section  C.l  we  show  that  they  have  the  following  functional  representations  [75] 

(2.11)  S*jk=e(6jk  +  {xj,Xk)i),  hjk  =  (AXj,Xk)i,  A  =  (-A)-1(<9r  -  u  •  V), 

where  (f,h) i  =  (V/  •  V/i)  is  a  Sobolev-type  sesquilinear  inner-product  [63]  and  the  operator  (—A)-1  is 
based  on  convolution  with  respect  to  the  Green’s  function  for  the  Laplacian  A  [90] .  Since  the  function  Xj  is 
real-valued  we  have  (XjiXk) l  =  (Xk,  Xj)i-,  which  implies  that  S*  is  a  symmetric  matrix.  The  function  A\j 
is  also  real- valued.  We  establish  in  Section  C.l  that  the  operator  A  is  skew-symmetric  on  a  suitable  Hilbert 
space,  which  implies  that  Akj  =  (Axk,Xj) l  =  ~(Xk,Axj) l  =  —{Axj,Xk)i  =  ^A*fc  which,  in  turn,  implies 
that  A*  is  an  antisymmetric  matrix,  hence  Akk  =  {Axk,Xk)i  =  0. 

Applying  the  linear  operator  (—A)-1  to  both  sides  of  the  cell  problem  in  equation  (2.8)  yields  the 
following  resolvent  formula  for  Xj 

(2.12)  Xj  =  (e  +  A)~1gj,  g3  =  (-A)"1^-. 


From  equations  (2.11)  and  (2.12)  we  have  the  following  functional  formulas  for  S*fc  and  A*k  involving  the 
antisymmetric  operator  A 


(2.13)  S*k  —  e  (djk  +  {(e  +  A)  1g3,(e  +  A)  1gk) i)  ,  A*k  —  {A{e  +  A)  1g^{e-\-A)  1gk)i- 

Since  A  is  a  skew-symmetric  operator,  it  can  be  written  as  A  =  iM  where  M  is  a  symmetric  operator.  We 
demonstrate  in  Section  C.l  that  M  is  self-adjoint  on  an  appropriate,  densely  defined  subset  of  a  Hilbert 
space. 

The  spectral  theorem  for  self-adjoint  operators  states  that  there  is  a  one-to-one  correspondence  between 
the  self-adjoint  operator  M  and  a  family  of  self-adjoint  projection  operators  {Q(A)}age  —  the  resolution 
of  the  identity  —  that  satisfies  [92]  limx-^mfE  Q(A)  =  0  and  limA_>SUp£  Q(X)  =  I  [92],  where  I  is  the 
identity  operator  and  E  is  the  spectrum  of  M.  Define  the  complex  valued  function  g,jk( A)  =  {Q(X)g3l  gk)i, 
j,  k  =  1, . . . ,  d,  where  gj  =  (—A )~1Uj  is  defined  in  (2.12).  Consider  the  positive  measure  g,kk  and  the  signed 
measures  Re njk  and  Ivngijk  associated  with  /j,jk( A),  introduced  in  equation  (A. 3).  Then,  given  certain 
regularity  conditions  on  the  components  u3  of  the  fluid  velocity  field  u ,  the  functional  formulas  for  5*k  and 
A*k  in  (2.13)  have  the  following  Radon-Stieltjes  integral  representations,  for  all  0  <  e  <  oo  (see  Section  C.l 
for  details) 


(2.14) 


Sh  =  £ 


Ijk 


dReg,jk(X)\ 
~  A2 


£2  +  . 


I  ) 


A  dim  g,jk(X) 
e2  +  A2 


A  key  feature  of  equation  (2.14)  is  that  parameter  information  in  e  is  separated  from  the  complicated 
geometry  and  dynamics  of  the  time-dependent  flow,  which  are  encoded  in  the  spectral  measure  g,jk.  This 
important  property  of  the  integrals  in  (2.14)  follows  from  the  non-dimensionalization  of  the  advection- 
diffusion  equation  discussed  in  the  paragraph  leading  to  equation  (2.2),  yielding  a  spectral  measure  p:jk 
that  is  independent  of  the  molecular  diffusivity  e.  An  alternate  formulation  of  the  effective  parameter 
problem  for  advection-diffusion  by  time-dependent  flows  was  discussed  in  [6],  which  used  a  different  non- 
dimensionalization,  yielding  a  Stieltjes  integral  representation  for  S£fc  involving  the  Peclet  number  Pe  of 
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the  flow  and  a  spectral  measure  that  depends  on  the  Strouhal  number.  However,  as  pointed  out  in  [18], 
the  Strouhal  number  dependence  of  the  measure  led  to  an  implicit  dependence  of  the  spectral  measure  on 
Pe.  This  restricts  the  utility  of  the  integral  representations,  such  as  rigorous  bounds  [7,  38]  which  depend 
explicitly  on  Pe  but  also  implicitly  on  Pe  through  the  moments  of  the  measure.  Our  formulation  has  no 
such  restrictions. 

The  periodic  homogenization  theorem  summarized  by  equations  (2.2)-(2.9),  as  well  as  its  many  varia¬ 
tions  [9,  72,  13,  15,  61,  5,  75,  76,  77,  62,  57],  depend  on  the  detailed  nature  of  the  fluid  velocity  field  u. 
They  also  depend  on  the  temporal  scaling  used  [15,  75,  57],  i.e. ,  what  value  of  7  is  used  in  equation  (2.5). 
However,  the  mathematical  structure  of  the  cell  problem  in  (2.8)  and  the  functional  form  of  D*  displayed  in 
equation  (2.9)  remain  unchanged  for  the  space-time  periodic  setting.  One  of  the  key  goals  of  the  present  work 
is  to  develop  a  rigorous  mathematical  framework  that  provides  the  Stieltjes  integral  representations  for  effec¬ 
tive  diffusivity  tensor  D*  displayed  in  (2.14),  for  space-time  periodic  u.  Due  to  the  time-dependence  of  the 
fluid  velocity  field,  one  must  employ  the  spectral  theory  of  unbounded  self-adjoint  operators,  which  is  much 
more  subtle  and  challenging  than  that  of  bounded  operators.  We  will  demonstrate  that  this  mathematical 
framework  depends  only  on  the  structure  of  the  cell  problem  in  (2.8)  and  the  presence  of  an  inner-product 
in  the  functional  form  of  D*  in  (2.9).  In  particular,  the  theoretical  development  is  insensitive  to  the  detailed 
nature  of  u  and  depends  only  on  its  boundedness  properties  (See  Corollary  C.2  for  details).  Consequently, 
our  results  given  here  apply  in  many  of  the  well  studied  systems  and  will  likely  apply  to  many  of  the 
homogenization  results  of  the  future. 

In  order  to  illustrate  the  rich  behaviors  that  can  arise  in  the  effective  diffusivity  tensor  D*  for  more 
general  velocity  fields  and  alternate  temporal  scalings,  we  now  briefly  discuss  some  key  variations  of  the 
theory  described  above.  When  the  fluid  velocity  field  is  mean-zero,  as  discussed  above,  then  equation  (2.7) 
holds  and  the  effective  diffusivity  tensor  D*  defined  in  (2.9)  is  constant  [57].  Consequently,  only  the  symmetric 
part  of  D*  plays  a  role  in  the  effective  transport  equation  displayed  in  (2.3).  Now  consider  a  more  general 
fluid  velocity  field 

(2.15)  u(t,  x)  =  Sau0(S1t,  Sx)  +  «i(t,  a:),  a  =  1,  7  =  2, 

which  is  the  superposition  of  a  weak ,  large-scale  mean  flow  6uo(62t,  Sx)  that  varies  on  large  spatial  and  slow 
time  scales,  with  a  mean-zero  periodic  flow  Ui(t,  x)  that  rapidly  fluctuates  in  space  and  time  [57].  If  Uo(t,  x) 
is  smooth  and  bounded,  the  homogenization  theorem  for  purely  periodic  velocity  fields  discussed  above  can 
be  rigorously  extended  to  the  present  setting  and  the  effective  transport  equation  in  (2.3)  is  replaced  by  [57] 

(2.16)  dt4>(t,  x)  =  u0(t,x)  •  Vcj)(t,  x)  +  V  •  [D*V(j)(t,  x)],  ^(0,  x)  =  (f>o(x), 

which  includes  an  advective  enhancement  in  transport  by  the  large-scale  mean  flow  Uq  [57].  In  this  case, 
the  effective  diffusivity  tensor  D*  is  completely  independent  of  the  mean  flow  u 0,  and  is  determined  by  the 
same  formula  in  equation  (2.9)  and  the  same  cell  problem  in  (2.8)  with  u  — >  U\  [57].  Consequently,  D*  is 
again  constant  and  only  the  symmetric  part  of  D*  plays  a  role  in  the  effective  transport  equation  displayed 
in  (2.16). 

This  problem  was  studied  in  [75]  for  scalings  in  (2.15)  different  than  a  =  1  and  7  =  2.  The  parameter 
a  determines  the  strength  of  the  mean  flow  tt0  relative  to  the  small  scale  periodic  fluctuations  U\.  When 
the  mean  flow  is  weak  compared  to  the  fluctuations,  to  leading  order,  D*  is  constant  and  independent  of  the 
mean  flow,  which  only  determines  the  transport  velocity  on  large  length  and  long  time  scales.  Consequently, 
only  the  symmetric  part  of  D*  plays  a  role  in  the  effective  transport  equation,  which  is  similar  to  that 
in  (2.16)  [75].  Regardless  of  the  values  of  a  and  7,  in  the  weak  mean  flow  regime,  the  components  D*fc  of 
the  effective  diffusivity  tensor  are  given  by  a  formula  analogous  to  equation  (2.9)  and  the  structure  of  the 
cell  problem  is  analogous  to  equation  (2.8).  There  are  three  distinct  behaviors  that  arise  as  the  values  of  a 
and  7  vary,  and  the  function  Xj  in  (2-8)  can  be  time-dependent  or  time-independent  ( dTXj  =  0)  [75]. 

As  we  discussed  in  Section  1,  the  constancy  of  the  effective  diffusivity  tensor  D*  is  not  consistent  with 
measurements  and  numerical  simulations  of  passive  tracer  transport  in  the  ocean  and  atmosphere.  However, 
when  the  fluid  velocity  is  active  on  both  the  slow  and  fast  time  scales,  with  7=1,  and  the  mean  flow  is 
equal  in  strength  or  stronger  than  the  periodic  fluctuations,  then  D*  is  a  function  of  space  and  time  [75]. 
Consequently,  in  the  effective  transport  equation,  the  antisymmetric  part  of  D*  contributes  to  a  purely 
rotational  (divergence-free)  enhancement  in  advective  transport,  while  the  symmetric  part  of  D*  contributes 
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to  an  enhancement  in  advective  and  diffusive  transport  [75].  This  is  consistent  with  observations  and  direct 
numerical  simulations  of  geophysical  flows  in  the  climate  system. 

In  Section  C.l  we  provide  a  mathematically  rigorous  framework  that  leads  to  the  Stieltjes  integral 
representations  in  (2.14),  for  both  the  symmetric  and  antisymmetric  parts  of  the  effective  diffusivity  tensor 
D*.  This  formulation  is  based  on  the  spectral  theorem  for  unbounded  self-adjoint  operators  in  Hilbert  space, 
which  is  based  on  an  axiomatic  construction  of  Hilbert  space.  Consequently,  the  integral  representations 
for  D*  depend  only  on  abstract  properties  of  the  underlying  self-adjoint  operator  and,  in  particular,  on 
boundedness  properties  shared  by  a  large  class  of  fluid  velocity  fields  u ,  including  all  those  discussed  in 
this  section.  In  Section  A,  we  review  the  spectral  theory  of  unbounded  operators.  In  Section  C  we  give 
two  natural  Hilbert  space  formulations  of  the  effective  parameter  problem  for  D*  which  lead  to  its  promised 
integral  representations.  In  Section  D  we  use  powerful  methods  of  functional  analysis  to  prove  that  the  two 
formulations  are  equivalent  and  discuss  the  theoretical  and  computational  advantages  of  each  approach. 

3.  Fourier  methods.  In  this  section  we  discuss  how  Fourier  methods  can  be  employed  to  convert  the 
eigenvalue  problem  Atpi  =  iX iipi,  A;  £  1, 1  £  N,  into  an  infinite  set  of  algebraic  equations  involving  the  Fourier 
coefficients  of  the  eigenfunction  ipi .  This  will  be  used  in  Section  4  to  compute  the  discrete  component  of  the 
spectral  measure  y,jk  underlying  the  integral  representations  for  the  effective  diffusivity  tensor  D*  displayed 
in  equations  (2.14)  and  (D.13).  For  notational  simplicity,  we  set  u  — >  u  and  (t, £)  (t,x)  and  write  the 

operator  A  in  (2.11)  as  A  =  (— A)-1(<9*  +  «•  V).  We  will  focus  on  the  fluid  velocity  field  u  in  equation  (1.1), 
although  the  methods  discussed  here  extend  to  a  large  class  of  fluid  velocity  fields,  namely,  those  expressible 
as  a  finite  sum  of  Fourier  modes. 

In  Section  C.l  we  showed  that  M  =  —iA  is  a  self-adjoint  operator  on  an  appropriate,  densely  defined 
subset  of  a  Hilbert  space  3^  C  L2(T  x  V)  given  in  equation  (C.3).  Since  the  space-time  periodic  velocity 
field  u  in  (1.1)  is  mean-zero,  we  require  that  the  elements  of  34?  are  also  mean-zero,  i.e. ,  /  £  34?  implies 
that  (/)  =  0.  In  Section  D.l  we  discuss  how  a  decomposition  of  the  Hilbert  space  34?  provides  a  natural 
decomposition  of  the  spectral  measure  fj,jk  into  its  discrete  and  continuous  components,  as  displayed  in 
equation  (D.8).  Here,  we  demonstrate  how  this  mathematical  framework  can  be  used  to  compute  the 
discrete  component  of  y jk ,  hence  the  discrete  component  of  the  Stieltjes  integral  representation  for  the 
effective  diffusivity  tensor  D*  displayed  in  equation  (D.13). 

Toward  this  goal,  rewrite  the  eigenvalue  problem  Aipi  =  iXiipi  as 

(3.1)  (dt  +u  ■  V)ipl  =  -iXiAip,,. 

For  spatial  dimension  d  =  2,  the  set  of  functions  {exp[j(£f  +  mx  +  ny)\  \  £,m,n  £  Z}  are  complete  in 
L2(T  x  V),  hence  the  Hilbert  space  34?.  Consequently,  the  eigenfunction  tpi  £  34?  can  be  represented  as 

(3.2)  <Pl(f,x)=  <m,n^U+mX+nV\  4,m,n  =  (w-e,(fl+”+nS))’ 

f,m,nGZ 


where  we  have  denoted  x  =  (x,y)  and  (•,•}  is  the  L2(T  x  V)-inner-product.  Write  the  fluid  velocity  field 
in  (1.1)  as  u(t,x )  =  ui(x)  +  dcost  U2(x),  where  u\  =  (cosy, cosir)  and  «2  =  (siny, since),  and  write 


(3.3)  cos t  sin y  =  -Uelt  +  e~lt)(eiy  -  e~iy)  =  -  el(t~y)  +  el(~t+y) 

4*  4 1 

Consequently,  substituting  the  formula  in  (3.2)  for  ipi  into  equation  (3.1)  yields 


-  el(-*-y)) 


(3.4) 


gi(£t-\-mx-\-ny) 


iXi(m2  +  n2)a\mn 


+,m((fl[m,n_1  +  a\,m  ,72+ 1  )/2 

+d(a{ 

—  1,772,72—1  ^—1,772,72+1  ^+1,772,72—1  ae+i  ,772,72+1  )/+)) 

+m((ae,m-l,n  +  aim+l,n)/2 

+0(a\_ 

l,m-l,n  a£-l,m+l,n  +  a£+l  ,772—1,72  a£+l,m+l,n)  /  (4*)) 
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Since  {el^t+T7l:E+ny^}  is  complete  in  the  underlying  Hilbert  space,  Yht  m  n  el('H+mx+nv^ci,m,n  =  0  only  if 
ce,m,n  =  0  for  all  t,m,n  £  Z  [92].  Consequently,  equation  (3.4)  yields  the  following  infinite  set  of  algebraic 
equations  for  the  Fourier  coefficients  a\  m  n  of  the  eigenfunction  ipi 

(3-5)  it  o\,m,n  +  2  {m(ae,m,n- 1  +  ai>,m,n+l)  +  n(at,m-l,n  +  a£,m+  l,n)] 

Q 

[^(^— l,m,n— 1  &£—  l,m,n+l)  ^{p,£— l,m— l,n  l,m+l,n) 

+™(aJ+l,m,n-l  -  ai+  l,m,n+l  )  +  n(al>+1 

,m—l,n  ^+l,m-t-l,n)] 

=  i\i(m2  +  n2)a\mn,  t,m,n£Z, 

The  solution  of  equation  (3.5)  determines  both  the  eigenvalues  A /  and  the  Fourier  coefficients  a\mn  of  the 
eigenfunction  ipi.  Notice  that  for  to  =  n  =  0,  equation  (3.5)  implies  that  it  a\  0  0  =  0  for  all  i  £  Z.  This  is 
consistent  with  the  requirement  that  ipi  is  mean-zero  in  space,  which  implies  that  a[  0  0  =0  for  all  t  £  Z. 
This  will  be  used  in  Section  4  to  compute  the  discrete  component  of  the  spectral  measure  pjk- 

Towards  this  goal,  we  now  determine  the  spectral  weights  ( gj ,  ipi) i  of  the  measure  pijk  displayed  in  (D.13) 
in  terms  of  the  Fourier  coefficients  aj  mn  of  ipi.  Recall  from  equation  (1.1)  that  u  =  (iq,  112)  with 


(3.6) 


u\ (t ,  x,  y)  =  cos  y  +  0  cos  t  sin  y 

=  I 


->  +  £( 


ei(t+y)  _  ej(-t+y)  _  gi(-t-y) 


)■ 


and  U2{t,x,y)  =  ui(t,y,x).  This  equation,  the  formula  =  (uj,ipi) 2  in  equation  (D.14),  the  Fourier 

expansion  of  ipi  in  (3.2),  and  the  orthogonality  of  the  set  { et(M+mx+ny)y,  imply  that 


(3.7) 


1  0 

{^£91)1  =  ^  (a0,0,l  +  a0,0,-l)  +  (al,0,l  —  al,0,-l  +  a— 1,0,1  —  a-l,0,-l) 

1  6 

(<7’2j52)i  =  ^  (a0,l,0  +  ao,-i,o)  +  ^  (al,l,0  —  °1,-1,0  +  1,1,0  —  °-l,-l,o) 


with  {gi,<pi)i  =  {<pi,gi)v 

When  6  =  0  in  the  velocity  field  of  equation  (1.1)  so  that  u  =  u(x)  is  time-independent,  the  operator 
A  no  longer  involves  the  time  derivative  and  the  associated  eigenfunction  ipi  =  ipi (x)  is  time- independent  as 
well.  In  this  case,  the  system  of  equations  in  (3.5)  reduces  to 


(3.8) 


TO 


—  (c 

O  V 


a. 


1 

m,n+ 1 


n ,  1 

)  +  7:(am-l,n  +  1 


'm+l,n 


)  =  A;  (to2  +  n2)c 


to,  n  £  Z, 


while  equation  (3.7)  for  the  reduces  to 

(3-9)  {<Pi,9i)i  =  \  (ao,i  +«o,-i)  >  (^2,32)1  =  ^  (ayo +  °-i,o)  ■ 

4.  Numerical  Results.  In  Section  3,  we  used  Fourier  methods  to  transform  the  eigenvalue  problem 
Aipi  =  iXitpi  for  the  operator  A  in  (2.11),  involving  the  fluid  velocity  field  in  (1.1),  into  an  infinite  system  of 
algebraic  equations  for  the  Fourier  coefficients  of  the  eigenfunctions  ipi.  In  this  section,  we  truncate  the  infinite 
system,  convert  it  to  a  generalized  eigenvalue  problem,  and  numerically  compute  the  discrete  component  of 
the  spectral  measure  underlying  the  integral  representations  for  the  symmetric  S*  and  anti-symmetric  A* 
parts  of  the  effective  diffusivity  tensor  D*,  displayed  in  equations  (2.14)  and  (D.13). 

By  restricting  the  indices,  —  M  <  £,m,  n  <  M ,  and  imposing  the  boundary  conditions 

(4.1)  aim,n  =  0  if  max(|4|TO|,|n|)  >M, 


the  infinite  systems  of  equations  in  (3.5)  and  (3.8)  become  finite  sets  of  equations.  Consider  the  fluid 
velocity  field  in  (1.1)  with  parameter  6  £  [0,1].  In  the  dynamic  ( 9  >  0)  and  steady  ( 9  =  0)  cases,  the 
bijective  mappings  Qd{t,m,n)  and  0s(m,  n)  defined  by 

(4.2)  Qd(t,  to,  n)  =  (M  +  to  +  1)  +  (M  +  n){2M  +  1)  +  (M  +  t)(2M  +  l)2, 

0s(m,  n)  —  ( M  +  to  +  1)  +  (M  +  n)(2M  +  1), 
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map  the  corresponding  finite  set  of  equations  to  matrix  equations.  In  particular,  they  become  generalized 
eigenvalue  problems 

(4.3)  B  ai  =  XiCai. 

Here  B  is  a  symmetric  matrix  and  C  is  a  diagonal  matrix  of  size  (2 M  +  l)3  x  (2 M  +  l)3  for  the  dynamic  case 
and  of  size  (2 M  +  l)2  x  (2 M  +  l)2  for  the  steady  case.  More  specifically,  B  is  Hermitian  in  the  dynamic  case 
and  is  real-symmetric  in  the  steady  case,  while  the  matrix  C  is  real-symmetric  and  diagonal  in  both  cases. 
Since  B  and  C  are  symmetric  matrices,  the  generalized  eigenvalues  A;  are  real- valued  and  the  eigen- vectors 
ai  -  consisting  of  the  Fourier  coefficients  for  ipi  -  satisfy  the  orthogonality  condition  [74] 

(4.4)  aJC  ak  =  6jk. 

The  matrix  C  is  positive  semidefinite  and  diagonal.  However,  in  the  steady  case,  in  both  of  the  matrices 
B  and  C,  the  row  and  column  associated  with  the  Fourier  coefficient  a}Q  0  is  all  zeros,  as  can  be  seen  from 
equation  (3.8).  Since  the  eigenfunction  c pi  is  mean-zero,  we  have  that  al0  0  =  0.  Therefore,  we  simply  remove 
this  row  and  column  in  both  B  and  C  so  that  C  becomes  strictly  positive  definite. 

In  the  dynamic  case,  the  entries  of  the  matrix  C  do  not  depend  on  £.  Consequently,  we  have  that 
C  =  diag(Cs, . . . ,  Cs),  where  we  have  denoted  by  Cs  the  matrix  C  in  the  steady  case.  Therefore,  there 
are  2 M  +  1  diagonal  entries  in  C  with  the  value  zero,  corresponding  to  m  =  n  =  0.  The  entries  of  the 
corresponding  rows  and  columns  of  the  matrix  B  are  all  zero  except  for  the  diagonal  entry,  which  has  the 
value  £,  as  can  be  seen  from  equation  (3.5).  This  implies  that  £aleo  0  =  0  for  all  —  M  <  £  <  M.  Since  the 
eigenfunction  ipi  is  mean-zero  in  time  and  space,  we  have  that  a\  0  0  =  0  for  all  —M  <  £  <  M,  which  is 
consistent  with  the  above  observation.  We  therefore  simply  remove  the  corresponding  rows  and  columns  in 
both  B  and  C  so  that  C  becomes  strictly  positive  definite.  This  method  of  removing  the  null  space  common 
to  both  B  and  C  is  called  deflation  [74]. 

Now  that  the  matrix  C  is  strictly  positive  definite  and  diagonal,  the  matrix  Cq  is  well  defined  for  all 
q  £  R,  with  entries  (C9)^-  =  C ^Sij,  where  is  the  ith  diagonal  of  the  matrix  C  raised  to  the  power  q , 
and  CqC~q  =  I.  Consequently,  the  generalized  eigenvalue  problem  in  equation  (4.3)  can  be  written  as  the 
following  standard  eigenvalue  problem 

(4.5)  C”1/2BC”1/2v;  =  A  ivi,  vi  =  C1/2at. 

Since  B  is  a  symmetric  matrix  and  C  is  diagonal,  the  matrix  C-1/2BC-1/2  is  also  symmetric  with  real¬ 
valued  eigenvalues  and  orthonormal  eigenvectors.  From  the  orthogonality  relation  =  Sjk  we  recover 

equation  (4.4)  via  Vi  =  C 1,/2a;  in  (4.5). 

In  summary,  our  numerical  method  is  the  following.  Create  the  matrices  B  and  C  according  to  equa¬ 
tion  (3.5)  or  (3.8)  and  the  corresponding  bijective  mapping  in  (4.2).  Remove  the  rows  and  columns  of  the 
matrices  B  and  C  corresponding  to  Qj  =  0.  Compute  the  eigenvalues  A;  and  eigenvectors  Vi  of  the  symmetric 
matrix  C~1//2BC_1/2.  The  computed  Fourier  coefficients  of  the  eigenfunction  tpi  are  given  by  a;  =  C_1/2u;. 
The  eigenvalues  associated  with  the  discrete  component  of  the  spectral  measure  displayed  in  equation  (D.13) 
are  given  by  A;,  while  the  spectral  measure  weights  (<?i,</?z)i  and  (g2,iPi)i  in  (D.13)  are  determined  by  the 
vector  ai  via  equation  (3.7)  or  (3.9). 

In  our  computations,  we  used  for  the  steady  case  M  =  150,  yielding  matrices  of  size  (2 M  +  l)2  —  1  = 
90,  600,  while  in  the  dynamic  case  we  used  M  =  20,  yielding  matrices  of  size  (2 M  + 1)3  —  (2 M  + 1)  =  68,  880. 
The  eigenvalues  and  eigenvectors  of  the  symmetric  matrix  C-1/2BC-1/2  were  computed  using  the  Matlab 
function  eiyf).  The  stability  of  the  computations  are  measured  in  terms  of  the  condition  numbers  ICi  of  the 
eigenvalues  A i,  which  are  the  reciprocals  of  the  cosines  of  the  angles  between  the  left  and  right  eigenvectors. 
Eigenvalue  condition  numbers  close  to  1  indicate  a  stable  computation.  Our  eigenvalue  computations  are 
extremely  stable  with  max;  |1  —  /C;|  ~  10-14,  which  were  computed  using  the  Matlab  function  condeig(). 

Displayed  in  Fig.  1  are  our  computations  of  the  discrete  component  of  the  spectral  measure  d/rn(A)  = 
(dA)  associated  with  the  fluid  velocity  field  u  displayed  in  equation  (1.1),  for  (a)  the  steady 
( 9  =  0)  and  (b)  the  dynamic  {9  =  1)  settings.  Here,  the  spectral  weights  =  \{gi,(pi)i\2  are  determined 

by  equations  (3.9)  and  (3.7),  respectively.  Consistent  with  the  symmetries  of  the  flows  [16],  we  have  Hu  =  ^22, 
while  Re  g,\2  =  0  and  Im  /112  =  0,  up  to  numerical  accuracy  and  finite  size  effects. 
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Fig.  1.  Computations  of  spectral  measures  and  effective  diffusivities  for  steady  and  dynamic  flows.  The  spectral  measure 
/in  associated  with  the  flow  in  (1.1)  are  displayed  for  (a)  the  steady  setting  and  (b)  the  dynamic  setting  with  the  associated 
effective  diffusivity  displayed  in  (c)  and  (d),  respectively.  In  the  steady  case  (a),  the  limit  point  of  the  measure  near  A  =  0 
has  small  measure  mass  with  mu  <  10— 30,  leading  to  the  asymptotic  behavior  SJj  ~  e1/2  for  £  <  1,  displayed  in  (c).  In  the 
dynamic  case  (b),  the  significant  measure  mass  mu  >  10“ 10  near  A  =  0  leads  to  the  asymptotic  behavior  SjJh  ~  1  for  g  <  1, 
displayed  in  (d). 


For  the  2D  steady  cell  flow  in  (1.1)  with  9  =  0,  it  is  known  [31,  71]  that  S],  ~  e1/2  for  e  <C  1.  Our 
computation  of  S]-,  displayed  in  Fig.  1(c)  is  in  excellent  agreement  with  this  result,  with  a  computed  critical 
exponent  of  ss  0.52  having  an  error  of  only  4%  relative  to  its  true  value  0.5.  Reducing  M  from  150  to 
100  changes  the  value  of  the  critical  exponent  by  less  than  0.0015,  indicating  that  the  value  of  M  =  150 
is  sufficiently  large.  In  this  steady  setting,  the  underlying  operator  (— A)_1[mi  •  V]  is  compact  [14]  and 
therefore  has  bounded,  discrete  spectrum  away  from  the  spectral  origin,  with  a  limit  point  at  A  =  0  [90]. 
The  limit  point  behavior  of  the  measure  /in  can  be  seen  in  the  rightmost  panel  of  Fig.  1(a).  The  decay 
of  for  vanishing  e  is  due  to  the  magnitude  of  the  measure  masses  mu(()  <  10-30  for  |Aj|  -C  1,  with  a 
significant  spectral  gap  near  the  limit  point.  The  rigorous  result  [31,  71]  ~  e1^2  as  e  — >  0  reveals  that 

the  spectrum  of  the  operator  (— A)-1[iti  •  V]  at  A  =  0  is  either  continuous  or  it  is  discrete  with  zero  mass, 
otherwise  S']-,  would  diverge  as  e  — >  0. 

In  contrast,  as  shown  in  Fig.  1(b),  the  spectral  measure  pn  associated  with  the  time-dependent  fluid 
velocity  field  in  (1.1),  with  9  =  1,  has  significant  values  of  mn(l)  near  the  spectral  origin,  with  toh(Z)  >  10_1° 
more  than  20  orders  of  magnitude  greater  than  that  of  the  steady  flow.  A  limit  point  behavior  in  the  measure 
pn  near  A  =  0  can  be  seen  in  the  rightmost  panel  of  Fig.  1(b).  It  is  interesting  to  note  that  the  support  of  the 
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measure  supp  pn  satisfies  supp  pn  C  [— M,M ]  for  all  values  of  M  investigated.  Due  to  the  significant  mass 
of  the  measure  near  the  spectral  origin  and  its  uniform  nature,  as  shown  in  the  center  panel  of  Fig.  1(b),  the 
effective  diffusivity  has  an  0(1)  behavior,  Sjd  ~  1  for  e  <C  1,  as  shown  in  Fig.  1(d).  This  is  consistent  with 
numerical  computations  of  Sjd  using  alternate  methods  [16].  This  0(1)  behavior  of  S*x  has  been  attributed 
to  the  Lagrangian  chaotic  behavior  of  the  flow  [16,  104]. 

Appendix  A.  Spectral  theory  of  unbounded  self-adjoint  operators  in  Hilbert  space.  The 

theory  of  unbounded  operators  in  Hilbert  space  was  developed  largely  by  John  von  Neumann  and  Marshall 
H.  Stone.  It  is  considerably  more  technical  and  challenging  than  that  of  bounded  operators,  as  unbounded 
operators  do  not  form  an  algebra,  nor  even  a  linear  space,  because  each  one  is  defined  on  its  own  domain. 
In  this  section,  we  review  the  spectral  theory  for  such  operators  and,  in  particular,  the  celebrated  spectral 
theorem  for  self-adjoint  operators  [84,  92]. 

An  operator  is  not  determined  unless  its  domain  is  known.  Let  <f>i  and  $2  be  operators  acting  on  a 
Hilbert  space  J4?  with  domains  D($i)  and  D($ 2),  respectively,  £>(<!>,)  C  J(?,  i  =  1,2.  They  are  said  to  be 
identical ,  in  symbols  $1  =  $2,  if  and  only  if  H($i)  =  D(4> 2)  and  $1/  =  $2/  for  every  /  of  their  common 
domain.  They  are  said  to  be  equal  in  the  set  ,  in  symbols  $1  =  $2,  if  and  only  if  X/  C  D(<f>i)  nD)^)  and 
•I>i/  =  $2/  for  every  /  £  X/ .  The  operator  $2  is  said  to  be  an  extension  ( proper  extension )  of  the  operator 
$1  if  D(<I>i)  C  H($2)  (H($i)  C  D(' F2))  and  the  operators  $2  and  $1  are  equal  in  D(4>i)  [92], 

Consider  the  sesquilinear  inner-product  (•,•}  associated  with  satisfying  (aip,bp)  =  ab(ip,p)  and 
(ip,  p)  =  (ip,  ip)  for  all  ip,p  £  Jl?  and  a,  b  £  C,  where  z  denotes  complex  conjugation  of  z  £  C.  The  .^-inner- 
product  induces  a  norm  ||  •  ||  defined  by  ||^||  =  (ip,  ip)1/2.  A  linear  operator  $  is  said  to  be  closed  if,  whenever 
the  sequences  {/„}  and  {$/„}  exist  and  converge  in  the  norm  ||  •  ||  to  the  limits  /  and  h,  then  /  £  D(4>)  and 
<h/  =  h  [92].  The  (Hilbert  space)  adjoint  <f>*  of  $  is  defined  by  ($>ip,  p)  =  (ip,  §*p)  for  every  ip  £  £>($)  and 
p  £  £)($*).  The  adjoint  <f>*  of  $  is  uniquely  determined  when  the  domain  D(<f>)  determines  Jtf,  i.e.,  the 
smallest  closed  linear  manifold  containing  D($)  is  the  Hilbert  space  Xtf?  [92].  In  this  case,  D($)  C  £)($*) 
and  <h*  is  a  closed  linear  operator  [92],  The  operator  <f>  is  said  to  be  symmetric  if  $  =  4>*.  The  operator  <f>  is 
said  to  be  self-adjoint  if  $  =  <f>*.  A  symmetric  operator  is  said  to  be  maximal  if  it  has  no  proper  symmetric 
extension.  A  self-adjoint  operator  is  a  maximal  symmetric  operator  [92]. 

The  operator  $  is  said  to  be  bounded  (in  operator  norm)  if  ||$||  =  sup^6^.  |ui|=1}  ||$V’II  <  00.  A 
bounded  linear  symmetric  operator  is  self-adjoint  if  and  only  if  its  domain  is  Jif  [92].  Conversely,  the 
Hcllinger-Toeplitz  theorem  states  that,  if  the  operator  $  satisfies  (&ip,  p)  =  (ip,  $p)  for  every  ip,p  £  Jtf, 
then  $  is  bounded  on  [84,  92].  This  indicates  that,  if  is  an  unbounded  symmetric  operator  on  Xrff,  then 
it  is  self-adjoint  only  on  a  densely  defined  proper  subset  of  Jtf. 

The  spectrum  E  of  a  self-adjoint  operator  $  on  a  Hilbert  space  Jtf  is  real- valued  [84,  92].  If  <f>  is  also 
bounded,  then  its  spectral  radius  equal  to  its  operator  norm  ||$||  [84],  i.e., 

(A.l)  SC[-||$||,||$||]. 

If  $  is  instead  unbounded,  its  spectrum  E  can  be  an  unbounded  subset  of,  or  can  even  coincide  with  the  set 
of  real  numbers  R  [92]. 

We  now  summarize  the  spectral  theorem  for  self-adjoint  operators  [92],  Let  $  be  a  self-adjoint  operator 
with  densely  defined  domain  D(<f>)  C  Jtf.  If  $  is  bounded  then  we  simply  take  D(4>)  =  Xtff.  The  spectral 
theorem  states  that  there  is  a  one-to-one  correspondence  between  the  self-adjoint  operator  $  and  a  family 
of  self-adjoint  projection  operators  {Q(A)}a£E  —  the  resolution  of  the  identity  —  that  satisfies  [92] 

(A.2)  lim  Q(A)  =  0,  lim  Q(A)  =  I, 

A — y  inf  S  A — ^  sup  S 

where  0  and  /  denote  the  null  and  identity  operators  on  Jf,  respectively.  Furthermore,  the  complex-valued 
function  of  the  spectral  variable  A  defined  by  p^v(X)  =  (Q(\)ip,p)  is  strictly  increasing  for  A  £  E  and  of 
bounded  variation  for  all  ip,  p  £  £)($)  [92]. 

By  the  sesquilinearity  of  the  inner-product  and  the  fact  that  the  projection  operator  <3(A)  is  self-adjoint, 
the  function  p^v(X)  satisfies  p^(X)  =  /Z^MA).  Moreover,  the  function  p^(X)  is  real- valued  and  positive 
MV’b(A)  =  (Q(X)ip,ip)  =  (Q(\)ip,Q(\)ip)  =  ||Q(A)V>  ||2  >  0.  Consider  the  associated  real- valued  functions 

R.e/q^v(A)  =  —  (pipip(X)  +  Abv(A))  ,  Im/i^,¥>(A)  =  —  (pj/j^( A)  —  /^^(A))  , 


(A.3) 


12 


N.  B.  Murphy  et  al. 


where  i  =  y/—l,  Re  p^(X)  =  p^(X)  and  Im/q^,(A)  =  0.  With  each  of  these  strictly  increasing  functions  of 
bounded  variation,  we  associate  Stieltjes  measures  [91,  92,  36] 

(A. 4)  dp^v(X)  =  d(Q(\)ip,(p),  dRep^v(X)  =  dRe  (Q(A)V’,  ¥>), 

dp^(X)  =  d||Q(A)V’  ||2-  dim n^v>(X)  =  dim  (Q (X)ip ,  ip) , 


which  we  will  denote  by  p^ip,  R ep$v,  and  Im/x^p.  We  stress  that  p^  is  a  positive  measure,  pj,v  is  a 
complex  measure,  while  Rep^  and  Im  p^v  are  signed  measures  [91,  92], 

The  spectral  theorem  also  provides  an  operational  calculus  in  Hilbert  space  which  yields  powerful  integral 
representations  involving  the  Stieltjes  measures  displayed  in  equation  (A. 4).  A  summary  of  the  relevant 
details  are  as  follows.  Let  F( A)  and  G( A)  be  arbitrary  complex-valued  functions  and  denote  by  St(F)  the 
set  of  all  ip  G  D(i E>)  such  that  F  €  L2(p^),  i.e. ,  F  is  square  integrable  on  the  set  £  with  respect  to  the 
positive  measure  p^,  and  similarly  define  @(G).  Then  @(F)  and  3>(G)  are  linear  manifolds  and  there  exist 
linear  operators  denoted  by  F(<F)  and  G(<&)  with  domains  @(F)  and  S>{G):  respectively,  which  are  defined 
in  terms  of  the  following  Radon-Stieltjes  integrals  [92] 


(A.5) 


(F(<i>)ip,ip) 

{F(*)i/>,G(*)<p) 


F{  A)  dp^v(X), 
F(X)G(X)  dp^(X), 


Vip  G  @{F),  ip  G  £($), 
Wipe  2>{F),  ip  G  f^(G), 


where  the  integration  in  (A.5)  is  over  the  spectrum  E  of  $  [84,  92]. 

The  mass  p^v  =  dp^(X)  of  the  Stieltjes  measure  p$v  satisfies  [92]  p^  =  limA_>sup E  p^v (A)  - 
lim^jnf  s  p^v(X).  Consequently,  equation  (A. 2)  and  the  Cauchy- Schwartz  inequality  yield 

/OO 

d(Q(X)ip,<p)  =  {ip, ip),  \plv\  <  Ill’ll  IMI  <  oo. 

-  OO 


Equation  (A. 6)  demonstrates  that  the  measures  in  (A. 4)  are  finite  measures ,  i.e.,  they  have  bounded  mass  [92]. 

Equation  (A.5)  can  be  generalized,  holding  with  suitable  notational  changes,  for  maximal  normal  oper¬ 
ators  [92].  Such  a  normal  operator  N  with  densely  defined  domain  D{ N)  C  commutes  with  its  adjoint 
N*,  i.e.,  NN*  =  N*N,  and  can  be  decomposed  as  N  =  $i  +  2,  where  $1  and  $2  are  self-adjoint  and 

commute.  The  spectrum  of  the  normal  operator  N  is  a  (possibly  unbounded)  subset  of  C  [92] .  A  special  case 
of  a  normal  operator  is  a  skew-adjoint  operator  satisfying  N*  =  — N.  It  can  be  decomposed  as  N  =  i<f>2  and 
since  $2  is  self-adjoint  having  purely  real  spectrum,  the  skew-adjoint  operator  N  =  f>2  has  purely  imagi¬ 
nary  spectrum  [92].  Consequently,  given  such  a  maximal  skew-adjoint  operator,  one  can  focus  attention  on 
the  self-adjoint  operator  $2  =  — ?N  without  having  to  resort  to  the  more  notationally  complicated  spectral 
theory  of  normal  operators. 

The  signed  measures  Rep^v  and  Imp^v  displayed  in  equation  (A. 4)  arise  naturally  when  considering  a 
maximal  skew-adjoint  operator  N  =  *<f>,  where  $  is  self-adjoint.  This  can  be  illustrated  by  considering  some 
special  cases.  Consider  the  functional  (F(ISl)ip,  G(N)</j)  involving  real-valued  Hilbert  space  members  F(N)ip 
and  G(N>,  so  that  (F(N)ip,  G(N)<p)  =  (G(N)<p,  F(N)ip)  G  M  and,  in  particular, 

(A. 7)  (F(N)^j,  G(N)ip)  =  i«F(NM  G(N)ip)  +  (G(N)y>,  F(N)ip)). 

Now  consider  the  special  cases  F{ N)  =  G(N)  and  F(N)  =  NG(N),  i.e.,  F(i\)  =  G(iX)  and  F(iX)  =  iXG(iX) 
in  equation  (A.5),  respectively.  It  follows  from  equations  (A.5)  and  (A. 7),  the  identities  Rex  =  (z  +  z)/2 
and  Imx  =  (z  —  z)/(2i),  and  the  linearity  properties  [92]  of  Stieltjes- Radon  integrals  with  respect  to  the 
functions  /q/,v(A)  and  jlpipi A)  that 

/OO 

|G(*A)|2dRe/x^(A), 

-OO 

/OO 

A  |G(*A)|2  dlm/q/^A). 

-00 
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An  important  property  of  a  self-adjoint  operator  <f>  which  will  be  used  later  is  that  its  domain  D(< f>) 
comprises  those  and  only  those  elements  ip  £  dff  such  that  the  Stieltjes  integral  A2  d/q^(A)  is  convergent. 
When  ip  £  _D(<I>)  the  element  is  determined  by  the  relations  [92] 

/oo  /*oo 

Ad/q/^A),  H^ll2  =  /  A2d/x^(A), 

-oo  J  — OO 


where  <p  is  an  arbitrary  element  in  D(<f>)  [92],  In  fact,  this  determines  the  one-to-one  correspondence  between 
the  self-adjoint  operator  $  and  its  resolution  of  the  identity  Q( A)  [92]. 

Appendix  B.  Time  derivative  as  a  maximal  normal  operator.  A  key  example  of  an  unbounded 
operator  is  the  time  derivative  dt  acting  on  the  space  L2  (T)  of  Lebesgue  measurable  functions  that  are  also 
square  integrablc  on  the  interval  T  =  [0,T],  say.  The  unboundedness  of  dt  as  an  operator  on  L2(T)  can  be 
understood  by  considering  the  orthonormal  set  of  functions  {<pn}  C  L2(T)  defined  by 


(B.l) 


<Pn(t)  =  P sin(mrt/T),  /3  =  y/2/T,  ( tpn ,  Wm)i  =  Snm, 


where  n,  m  £  N  and  (•,  -}2  denotes  the  sesquilinear  L2(T)-inner-product.  From  dt<pn  =  ( mrfi/T )  cos(nnt/T) 
and  ||9t</jn||2  =  ( nir/T )2,  it  follows  that  the  norm  of  the  members  of  the  set  {dt<pn}  grows  arbitrarily  large 
as  n  — >  oo.  This  clearly  demonstrates  the  unboundedness  of  the  operator  dt  with  domain  L2(T). 

When  one  also  imposes  periodic  or  Dirichlet  boundary  conditions,  simple  integration  by  parts  demon¬ 
strates  that  the  operator  dt  is  skew- symmetric  on  L2(jT)  so  that  —  idt  is  symmetric  with  respect  to  the 
sesquilinear  inner-product  (•,  -)2-  We  now  identify  an  everywhere  dense  subset  of  L2(T)  on  which  —idt  is  a 
bounded  linear  self-adjoint  operator  [84,  92].  Consider  the  class  of  all  functions  ip  £  L2(fT)  such  that  ip(t) 
is  absolutely  continuous  [85]  on  the  interval  T  and  has  a  derivative  ip' it)  belonging  to  L2(T),  i.e. ,  [92,  85] 


(B.2) 


=  \  *P  e  L2{T)  ip(t)  =c+  g(s)ds,  g  £  L2(T) 


where  the  constant  c  and  function  g(s)  are  arbitrary.  Now,  consider  the  set  .fi/7-  of  all  functions  tp  £  &/j- 
that  satisfy  the  periodic  boundary  condition  ^(0)  =  i.e.  functions  ^  satisfying  the  properties  of 

equation  (B.2)  with  g(s)ds  =  0.  In  order  to  help  clarify  the  ideas  that  were  discussed  in  Section  A  in 
terms  of  an  abstract  Hilbert  space  Jff,  we  also  consider  the  set  ^7-  of  all  functions  ip  6  that  satisfy  the 
Dirichlet  boundary  condition  ^>(0)  =  ip(T)  =  0,  i.e.  functions  ip  satisfying  the  properties  of  equation  (B.2) 
with  c  =  0  and  g(s)ds  =  0.  More  concisely, 

(B.3)  =  {ip  6  I  V’(O)  =  V’CO})  =  {ip  €  I  V'(O)  =  =  0}- 


These  function  spaces  satisfy  .2/7-  C  sti 7-  C  si 7-  and  are  each  everywhere  dense  in  L2(T)  [92].  Let  the 
operators  B,  B,  and  B  be  identified  as  —  idt  with  domains  si 7-,  si 7-,  and  si r,  respectively.  Then,  B  is  a 
closed  linear  symmetric  operator  with  the  adjoint  B*  =  B ,  and  the  operator  B  is  a  self-adjoint  extension  of 
B  [92].  In  symbols,  this  means  that  B  =  B*  on  s/j-  and  D{B)  =  D(B*)  =  stfq-,  i.e.,  B  =  B*  on  stf-/-.  This 
establishes  that  the  operator  — idt  with  domain  si 7-  is  self-adjoint,  hence  dt  is  a  maximal  skew-symmetric 
(normal)  operator  on  stf- j-.  The  operator  idt  on  si 7-  has  a  simple  point  spectrum,  consisting  of  eigenvalues 
A  =  2mr/T,  n  £  Z,  with  corresponding  eigenfunctions  exp(i2nirt/T)  [92], 

Appendix  C.  Hilbert  spaces,  resolvents,  and  integral  representations  of  the  effective  diffu- 
sivity.  In  this  section  we  formulate  a  spectral  theory  of  effective  diffusivities  for  space-time  periodic  flows. 
In  Section  C.l  we  address  an  approach  suggested  in  [75],  while  in  Section  C.2  we  address  an  approach  sug¬ 
gested  in  [6].  In  each  case,  we  provide  a  rigorous  mathematical  framework  which  leads  to  Stieltjes  integral 
representations  for  both  the  symmetric  S*  and  antisymmetric  A*  parts  of  the  effective  diffusivity  tensor  D*  for 
space-time  periodic  flows,  involving  a  spectral  measure  of  an  unbounded  self-adjoint  operator.  In  Section  D 
we  use  the  one-to-one  correspondence  between  a  self-adjoint  operator  and  its  resolution  of  the  identity  [92], 
discussed  in  the  paragraph  containing  equation  (A. 9),  to  establish  that  the  two  approaches  are  equivalent. 
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C.l.  Scalar  fields  and  the  effective  diffusivity.  In  this  section  we  provide  an  abstract  Hilbert  space 
formulation  of  the  effective  parameter  problem  for  advection  enhanced  diffusion  by  a  space-time  periodic  fluid 
velocity  field  Consider  the  following  sets  T  =  [0,  T]  and  V  =  x^=1[0,£]  which  define  the  space-time 

period  cell  T  x  V  for  u(t,x).  Now  consider  the  Hilbert  spaces  A2(T)  and  A2(V)  of  Lebesgue  measurable 
functions  over  the  complex  field  C  that  are  also  square  integrable  on  T  and  V,  respectively.  Define  the 
associated  Hilbert  spaces  My  and  M^, 

(C.l)  My  =  {if  &L2(T)\ip(t)  =  ip(t  +  T)},  Jfv  =  {^ei2(V)|^)  =  i)(x  +  fe3)}, 

for  all  j  =  1 .... ;  d,  where  the  are  standard  basis  vectors.  Denote  by  (•)  space-time  averaging  over  T  x  V. 
Now  define  the  Hilbert  space  Myy  =  My®  My  with  sesquilinear  inner-product  (•,  •)  given  by  (if),  ip)  =  (ipTp), 
with  (tp,ip)  =  {4’i‘p)-  The  ^7-y-inner-product  induces  a  norm  ||  •  ||  given  by  ||^>||  =  (ip, ip)1/2  [36]. 

In  equation  (B.3)  we  defined  the  space  si fy  of  absolutely  continuous  T-periodic  functions  with  time 
derivatives  belonging  to  My.  We  also  discussed  that  siy  is  an  everywhere  dense  subset  of  the  Hilbert  space 
My-  [92].  We  now  define  the  Sobolev  space  My ,  which  is  also  a  Hilbert  space  [14,  35,  63], 

(C.2)  =  {iP&Mv\(\ViP\'l)v  <00}, 

where  (-)y  denotes  spatial  averaging  over  V.  Finally,  consider  the  Hilbert  space  M  and  its  everywhere  dense 
subset  &  defined  by 

(C.3)  M’  =  ;Mr®Mf),  &  =  s/r®M$. 

Recalling  that  ip  •  ip  =  ip^tp,  the  sesquilinear  ^-inner-product  is  given  by  (ip,<p)  1  =  (V-0  •  V<p)  with 
associated  norm  ||  •  ||i  given  by  HV’lli  =  (|V'0|2)1'2.  We  stress  that  ip  €  &  implies  ||<9t?/’||i  <  00  and 
HV’Hi  <  00.  In  the  case  of  a  time- independent  fluid  velocity  field  u(x)  we  set  M  =  &  =  My . 

We  now  use  properties  of  the  Hilbert  space  M  to  obtain  functional  formulas  for  the  symmetric  S*  and 
antisymmetric  A*  parts  of  the  effective  diffusivity  tensor  D*  defined  in  equations  (2.9)  and  (2.10),  involving 
the  solution  \j  °f  the  cell  problem  in  equation  (2.8)  and  a  maximal  skew-symmetric  operator  A  on  & .  We 
then  derive  from  the  cell  problem  a  resolvent  formula  for  Xj  involving  the  operator  A.  The  spectral  theorem 
discussed  in  Section  A  then  yields  the  Stieltjes  integral  representations  for  S*  and  A*.  We  will  henceforth 
assume  that  Uj,  Xj  £  &  for  all  j  =  1, . . . ,  d.  We  stress  that,  since  in  each  term  of  equation  (2.8)  the  function 
Xj  is  acted  upon  by  a  differential  operator,  we  may  assume  without  loss  of  generality  that  (xj)  =  0. 
Applying  the  linear  operator  (—A)-1  to  both  sides  of  the  cell  problem  in  equation  (2.8)  yields 

(C.4)  (-A)~1Uj  =  (e  +  A)xj, 

where  we  have  defined  A  =  (— A)_1(di  —  u- V).  The  operator  (—A)-1  is  based  on  convolution  with  respect  to 
the  Green’s  function  for  the  Laplacian  A  and  is  bounded  on  L2(V)  [90],  hence  My  .  Now  write  the  functional 
(ujXk)  hi  equation  (2.9)  as  [75] 

(C.5)  (uj Xk)  =  ([AA_1Uj]  Xk)  =  -(VA_1Uj  •  Vyfc)  =  ((-A)~1Uj,  Xk) i- 

This  calculation  will  be  rigorously  justified  in  Theorem  C.l  below.  Substituting  the  formula  for  (— A)_1itj 
in  (C.4)  into  equation  (C.5)  yields  equation  (2.11),  which  provides  functional  formulas  for  the  components 
S*fc  and  A*fc,  j,k  =  1  ,...,d,  of  S*  and  A*.  Equation  (C.4)  leads  to  the  resolvent  formula  displayed  in 
equation  (2.12).  From  equations  (2.11)  and  (2.12)  we  have  the  functional  formulas  for  S*k  and  A*fc  displayed 
in  equation  (2.13),  involving  the  operator  A.  The  following  theorem  establishes  the  promised  Stieltjes  integral 
representations  for  the  functional  formulas  for  S*fc  and  A*k  in  (2.13). 

Theorem  C.l.  The  operator  A  =  (— A)-1(<9t  —  u  •  V)  displayed  in  equation  (2.11)  is  a  maximal 
(skew- symmetric)  normal  operator  on  the  function  space  &  defined  in  equation  (C.3),  hence  M  =  —1A  is 
self-adjoint  on  fp .  Let  Q( A)  be  the  resolution  of  the  identity  in  one-to-one  correspondence  with  M .  Define  the 
complex  valued  function  Pjk(A)  =  (Q(A)gj ,  gk)\,  j,k  =  1, . . . ,  d,  where  gj  =  (— A)~1Uj,  as  defined  in  (2.12), 
and  (-,-)i  is  the  M -inner-product.  Consider  the  positive  measure  pkk  o,nd  the  signed  measures  Repjk  and 
Impjk  associated  with  g.jk( A),  introduced  in  equation  (A. 3).  Then,  for  Uj,Xj  £  &  and  all  0  <  e  <  00,  the 
functional  formulas  for  S*fc  and  A*fc  in  (2.13)  have  the  Radon-Stieltjes  integral  representations  displayed  in 
equation  (2.14). 


Spectral  analysis  of  advection-diffusion  by  dynamic  flows 
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Proof  of  Theorem  C.l.  We  first  establish  that  M  =  —  iA  is  a  self-adjoint  operator  on  JC  The  Sobolev 
space  in  (C.2)  is  the  closure  in  the  norm  (|Vi/’|2)v  of  the  space  of  all  twice  continuously  differentiable 
periodic  functions  in  j£y ,  and  all  the  elements  of  are  those  elements  of  j£y  which  have  square  integrable 
gradients  on  the  set  V  [14].  Furthermore,  the  elements  of  s/f  are  those  elements  of  Mj-  that  are  differentiable 
almost  everywhere  (except  on  a  set  of  Lebesgue  measure  zero),  have  square  integrable  derivatives  on  the 
interval  T ,  and  are  indefinite  integrals  of  their  derivative,  hence  continuous  [85].  Consequently,  /  £  & 
implies  that  [92,  85] 

(C.6)  ll/lloo  =  sup  |/(i,  a;) |  <  oo, 

(t,®)erxv 

almost  everywhere.  For  Uj  £  &  and  fixed  t  £  T,  equation  (C.6)  implies  that  [u(t,  ■)  •  V]  :  Jtfy  — >•  Jtfy, 
while  (—A)-1  :  [14].  In  particular,  for  /,  h  £  we  have  that  ((— A)-1/,  K)\  —  (f,h)  [14].  This 

justifies  the  calculation  in  equation  (C.5)  (see  also  the  discussion  leading  up  to  equation  (D.2). 

We  have  already  established  in  Section  B  that  the  operator  —  idt  with  domain  s/j-  is  self-adjoint  [92]. 
The  integral  operator  (—A)-1  is  self-adjoint  and  compact  on  [90].  Since  they  commute  on  s/j-^J^  [36], 
it  follows  that  the  operator  —i(—A )~ldt  is  self-adjoint  with  domain  s/q -®^£y,  hence  (— A)_19t  is  a  maximal 
(skew-symmetric)  normal  operator  on  the  same  domain  [92] . 

We  now  establish  that  the  operator  (— A)_1[it- V]  is  antisymmetric  and  compact  on  The  antisymme¬ 
try  of  this  operator  depends  on  the  incompressibility,  V  -u  =  0,  of  the  fluid  velocity  field  and  was  established 
in  [14,  75].  Since  the  operator  (—A)-1  is  compact  on  j£y  [90],  we  need  only  show  that  the  operator  u  ■  V 
is  bounded  on  This  is  established  by  the  following  calculation.  For  Uj,  f  £  & ,  equation  (C.6)  yields 

(C.7)  ||tx-V/||2  =  |(u.V/,u-V/)| 

—  \(ujdjfiukdkf) |  (triangle  inequality) 

jk 

<  m&X  Halloo  \(®jfi®kf)\ 

3  jk 

<  max||uj||^5Z  \\djf\\  ||3fe/||  (Cauchy-Schwartz) 

3  jk 

=  ma-x  II II So  [  II ^ / II 

3 

<  d  max  ||«X  ||cb/||2  (Cauchy-Schwartz) 

j 

=  d  max  llujU^;,  ||/||?. 

3 

This  demonstrates  that  the  operator  norm  ||it  •  V||  has  the  upper  bound  ||it  •  V||  <  \[d  max.,-  ||Mjj|oo  <  00 
and  establishes  that  (—A)-1  [it  •  V]  is  a  compact  operator  on  & .  Since  (—A)-1  [u  •  V]  is  antisymmetric  and 
bounded  on  it  is  a  maximal  (skew-adjoint)  normal  operator  on  hence  the  operator  — 1{— A)-1  [it  •  V] 
is  self-adjoint  on  &  [92]. 

Denote  M  =  — iA ,  where  A  =  (— A)_1(9t  —  it  •  V).  Since  —  z(— A)_1[u  •  V]  is  self-adjoint  on  &  and 
— 1{— A)_19t  is  self-adjoint  on  si- 7-  ®  ^y ,  the  operator  M  is  self-adjoint  with  domain  D(M)  Z>  (s/j-  <g>  J^y)  fl 
&  =  &  [92], 

The  complex- valued  functions  involved  in  the  functional  formulas  for  S*fc  and  A*k  in  equation  (2.13)  are 
F( A)  =  (e+zA)-1  and  G(A)  =  zA(e+zA)_1.  For  all  0  <  e  <  00,  we  have  |T(A)|2  =  (e2  +  A2)-1  <  e~2  <  00  and 
|G(A)|2  =  A2(e2  +  A2)-1  <  1.  Since  Hkk  is  a  finite  measure  for  all  k  =  1, . . . ,  d,  as  shown  in  equation  (A. 6), 
we  therefore  have  that  /  £  3>(F)  and  /  £  S>(G)  for  all  /  £  D{M )  when  0  <  e  <  00.  Since  Uj  £  JF  and 
(—A)-1  is  a  bounded  operator  on  JF,  we  have  that  gj  =  (—A)~1Uj  £  JF.  We  note  that  in  the  mean-zero 
setting,  (uj)  =  0  and  the  Fubini-Tonelli  theorem  [36]  imply  that  we  also  have  (gj)  =  0.  The  conditions  of 
the  spectral  theorem  are  thus  satisfied.  Consequently,  the  integral  representations  in  equation  (A. 5)  hold  for 
the  functions  F( A)  and  G(A)  defined  above,  involving  the  complex  measure  fijk-  The  discussion  leading  to 
equation  (A. 8)  then  establishes  the  integral  representations  for  5*k  and  A*k  displayed  in  equation  (2.14). 
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It  is  worth  noting  that  from  equations  (A. 6)  and  (C.5),  the  mass  fj^k  of  the  measure  p,jk  is  given  by 
jjL °jk  =  (gj,gk)  1  =  {(— A)~1Uj,  uk)2l  where  (-,-)2  denotes  the  L2{T  x  V)  inner-product.  Since  (—A)-1  is  a 
self-adjoint  operator  on  L2(V)  [90],  the  spectral  theorem  demonstrates  that 

(C.8)  g°jk  =  ((-A  )~1uj,uk)2  =  J \d{Q(\)uj,uk)2- 

In  other  words,  the  mass  of  the  measure  p,jk  is  the  first  moment  of  the  spectral  measure  d(Q(X)uj,uk)2  of 
the  self-adjoint  operator  (—A)-1,  where  Q( A)  is  the  resolution  of  the  identity  in  one-to-one  correspondence 
with  (—A)-1.  This  completes  the  proof  of  Theorem  C.l  □. 

We  conclude  this  section  with  a  discussion  regarding  an  extension  of  Theorem  C.l  to  a  broader  class  of 
fluid  velocity  fields,  summarized  by  the  following  corollary. 

Corollary  C.2.  Theorem,  C.l  can  be  extended  to  the  following  class  ^  of  fluid  velocity  fields  u, 
having  components  Uj,  j  =  1 , ...  ,d, 

(C.9)  ^  =  {uj  £  srf'q-  ®  J^y  |  3  0  <  C  <  oo  such  that  ||  ( — A)-1  [it,  •  V]||  <  C}. 

Corollary  C.2  states  that  the  requirement  Uj  £  ^  =  s#x  ®  can  be  weakened  to  Uj  £  &?x  ®  such 
that  the  operator  (—A)-1  [it  •  V]  is  bounded  on  & .  The  set  ^  in  (C.9)  is  non-empty.  We  established  this 
in  (C.7),  showing  that  {uj  £  s#x  ®  |  ||Mjj|oo  <  oo}  C  as  (—A)-1  is  a  bounded  operator  on  j£y  [90]. 

This  extension  of  Theorem  C.l  allows  for  spatially  unbounded  flows  with  square  integrable  singularities. 

Proof  of  Corollary  C.2.  There  are  three  places  in  the  proof  of  Theorem  C.l  which  require  a  certain 
amount  of  regularity  in  the  components  Uj  of  the  fluid  velocity  field  it.  One  requirement  was  that  the 
operator  (—A)”1  [it  •  V]  be  bounded  on  &  so  that  A  is  a  maximal  (skew-symmetric)  normal  operator  on 
& .  Another  regularity  requirement  of  Uj  appeared  in  the  calculation  in  equation  (C.5).  The  functional 
(■ UjXk )  in  equation  (C.5)  is  well  defined  for  Uj,Xk  G  ©>  as  the  Cauchy- Schwartz  inequality  yields 
\{ujXk)  |  <  llMjllllXfcll  <  °o,  while  the  functional  ((— A)-1Uj,  Xk)i  in  (C.5)  is  well  defined  for  Uj  £  Jiff® 
and  Xk  €  Jtfj-  ®  as  (—A)-1  :  [14] .  However,  the  intermediate  step  ( UjXk )  —  ([AA_1iij]  Xk) 

required  that  A-1^  has  square  integrable  spatial  derivatives  of  order  two,  i.e. ,  Uj  £  Mf  ®  .  Although, 

after  the  integration  by  parts,  this  requirement  was  weakened  to  Uj  £  Jiff  ®  J£y.  The  final  regularity 
requirement  on  u7  was  in  the  conditions  of  the  spectral  theorem  in  (A. 5).  Namely,  that  (— A)~1Uj  £  D{A)  Z> 
TP,  as  well  as  (—A )~1Uj  £  S>(F)  and  (—A )~xUj  £  @(G)  for  F( A)  =  (e  +  zA)_1  and  G(A)  =  iX(e  +  *A)_1. 
However,  we  demonstrated  that  /  £  @{F)  and  /  £  S>(G)  for  all  f  £  XF .  Since  (—A)-1  :  -£  [14], 

we  only  require  that  Uj  £  s^x  ®  =x£y.  This  allows  for  spatially  unbounded  flows  with  square  integrable 
singularities.  An  exposition  of  the  specific  details  is  beyond  the  scope  of  the  current  work.  This  concludes 
our  proof  of  Corollary  C.2  □. 

C.2.  Curl-free  vector  fields  and  effective  diffusivity.  In  this  section  we  provide  a  rigorous  mathe¬ 
matical  framework  for  an  alternate  formulation  [6]  of  the  effective  parameter  problem  for  advection  enhanced 
diffusion  by  space-time  periodic  fluid  velocity  fields.  This  approach  provides  analogous  formulas  to  those 
displayed  in  equations  (2.11)-(2.14)  involving  the  curl-free  vector  field  Vy7  displayed  in  equation  (2.8)  and 
a  maximal  (skew-symmetric)  normal  operator  acting  on  a  suitable  Hilbert  space.  Towards  this  goal,  recall 
the  Hilbert  spaces  M f  and  J^y  given  in  equation  (C.l)  and  the  function  space  s^x  given  in  equation  (B.3). 
Now  define  their  d-dimensional  analogues  over  the  complex  field  C, 

(C.10)  TLx  =  Tty  =  O^JTy,  Ax  =  <s>j— \£^X- 

By  the  Helmholtz  theorem  [56,  12],  the  Hilbert  space  Tty  can  be  decomposed  into  mutually  orthogonal 
subspaces  of  curl- free  Ttx,  divergence- free  Tt.,  and  constant  Tto  vector  Helds,  with  tty  =  Ttx©Tt.  ©Tto- 
The  orthogonal  projectors  associated  with  this  decomposition  are  given  by  Tx  =  — V(— A)_1V* ,  T,  =  V  x 
(— A)-1  Vx  ,  and  To  =  (•),  respectively,  satisfying  I  =  Tx  +  T.  +  To  [31,  71,  66].  Here,  A  =  diag(A, . . . ,  A) 
is  the  vector  Laplacian  with  inverse  A-1  =  diag(A_1, . . . ,  A-1),  (•)  denotes  space-time  averaging  over  the 
period  cell  Tx  V,  and  I  is  the  identity  operator  on  Tty.  Due  to  the  curl-free  vector  field  at  the  heart 
of  the  cell  problem  in  equation  (2.8),  we  will  find  particular  use  of  the  Hilbert  space  7tx ,  which  we  define  as 

Ttx  =  {if  £  Tty  |  Tif  =  i/>  weakly},  T  =  —  V(— A)_1V- , 


(C.ll) 
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where  we  have  denoted  Tx  by  T  for  notational  simplicity.  Since  (—A)-1  is  self-adjoint  on  Hx  [90],  it  is  clear 
from  integration  by  parts  that  T  is  a  symmetric  operator  on  Hx ,  and  since  it  is  also  a  projection  operator,  it 
is  bounded  with  operator  norm  ||r||  =  1.  Thus  T  is  self-adjoint  on  Hx  [92,  84].  Analogous  to  equation  (C.3), 
we  define  the  Hilbert  space  H  and  its  everywhere  dense  subset  F , 

(C.12)  H  =  jtT  0  Hx,  F  =  Ar®Hx- 

Denote  by  ||  •  ||  the  norm  induced  by  the  the  sesquilinear  inner-product  (•,  •)  associated  with  the  Hilbert  space 
H,  defined  by  (if,<p)  =  (if  •  p)  with  (if,p)  =  (ip, if).  Here,  if  ■  ip  =  ifTp,  ifT  denotes  transposition  of 
the  vector  if,  and  Tp  denotes  component-wise  complex  conjugation,  with  if  •  if  =  \if\2 .  We  will  henceforth 
assume  that  u,  'Vxj  £  F.  In  the  case  of  a  steady  fluid  velocity  field  u(x ),  we  set  H  =  F  =  Hx  ■ 

Since  the  fluid  velocity  field  u  is  incompressible,  there  is  a  real  skew-symmetric  matrix  H(f,a;)  satisfy¬ 
ing  [4,  5] 

(C.13)  u  =  V  •  H,  Ht  =  -H, 

where  H T  denotes  transposition  of  the  matrix  H .  Since  u  e  F,  we  have  that  the  space-time  periodic 
components  Y\jk  of  the  matrix  H  have  square  integrable  spatial  derivatives  of  order  two  on  the  set  V  and 
are  absolutely  continuous  with  square  integrable  temporal  derivatives  on  the  set  T  [14].  Consequently,  as 
in  (C.6),  || Hjfc || oo  <  oo  hence  H  is  bounded  in  operator  norm  ||H||  <  oo  on  Hr®H\>-  Due  to  the  skew- 
symmetry  of  H,  we  have  the  identity  [V  •  H  ]  •  V/  =  V  •  [H  V/].  Using  this  identity  and  the  representation 
of  the  velocity  field  u  in  (C.13),  the  advection-diffusion  equation  in  (2.1)  can  be  written  as  a  diffusion 
equation  [31,  71], 

(C.14)  dt(f  =  V-DV(f,  cf(0,  x)  =  4>o(x),  D  =  ei  +  H, 

where  D (t,  x)  =  el  +  H(t,  x)  can  be  viewed  as  a  local  diffusivity  tensor  with  coefficients 
(C.15)  D  jk  =  eSjk  +  H  jk,  j,k  =  l,...,d. 

The  cell  problem  in  (2.8)  can  also  be  written  as  the  following  diffusion  equation  [31,  71] 

(C.16)  drXj  =  V?.[D(VjXj-+e,.)],  <vafc)  =  0,  D  =  el  +  H, 

where  (Vt%fc)  =  0  follows  from  the  periodicity  of  \k ■  We  stress  that  equation  (C.14)  involves  the  slow  ( t ,  x) 
and  fast  variables  (r,  £),  while  equation  (C.16)  involves  only  the  fast  variables.  For  notational  simplicity,  we 
will  drop  the  subscripts  £  displayed  in  equation  (C.16). 

We  now  recast  the  first  formula  in  equation  (C.16)  in  a  more  suggestive,  divergence  form.  Define  the 
operator  T  :  Aq — >  Hr  by  (T if)j  =  dTifj,  j  =  1, . . . ,  d.  For  /  €  &  we  have  [31,  71,  36,  35] 

(C.17)  V(A~1)dTf=  A^TV/, 

so  that  [31,  71]  dTXk  =  AA ~ldTXk  =  V  •  (A_1T)Vxfc-  Define  the  vector  field  Ek  =  Vyj,  +  ek  and  the 
operator  er  =  el  +  S,  where  S  =  (— A)_1T  +  H  and  in  the  case  of  a  steady  fluid  velocity  field  u(x)  we 
have  S  =  H  and  er  =  D.  With  these  definitions,  the  cell  problem  in  (C.16)  can  be  written  as  V  •  crEk  =  0, 
(Ek)  =  ek,  which  is  equivalent  to 

(C.18)  V  •  =  0,  VxEk  =  0,  Jk  =  crEk ,  (Ek)  =  ek,  cr  =  el  +  S. 

The  formulas  in  (C.18)  are  the  quasi-static  limit  of  Maxwell’s  equations  for  a  conductive  medium  [38,  66], 
where  Ek  and  Jk  are  the  local  electric  field  and  current  density,  respectively,  and  cr  is  the  local  conduc¬ 
tivity  tensor  of  the  medium.  In  the  analytic  continuation  method  for  composites  [38,  65,  10],  the  effective 
conductivity  tensor  cr*  is  defined  as 

(C.19)  (Jk)=cr*(Ek). 

The  linear  constitutive  relation  Jk  =  aEk  in  (C.18)  relates  the  local  intensity  and  flux,  while  that  in  (C.19) 
relates  the  mean  intensity  and  flux.  Due  to  the  skew-symmetry  of  S,  the  intensity-flux  relationship  in  (C.18) 
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is  similar  to  that  of  a  Hall  medium  [45,  31,  71].  The  precise  relationship  between  the  effective  parameters 
D*  and  er*  are  established  in  Lemma  C.4  below. 

Analogous  to  equation  (2.11),  the  components  S*fc  and  A*k,  j,k  =  of  the  symmetric  S*  and 

antisymmetric  A*  parts  of  the  effective  diffusivity  tensor  D*  can  be  represented  by  the  following  functional 
formulas  in  terms  of  the  T~L  inner-product  (•,•), 

(C.20)  S*k  =  e{Sjk  +  (Vxj,  Vxa,}),  A*jk  =  (AVXj,  VXfc),  A  =  TST. 

Equation  (C.20)  follows  from  the  formula  D*fc  =  eSjk  +  (uj  Xk)  in  equation  (2.9)  and  the  cell  problem  in 
equation  (C.18)  written  as  V  •  crVxj  =  — V  •  He.,  =  — Uj ,  yielding 

(C.21)  { uj  Xk)  =  -([V  •  crVxj], Xk) 

=  (o’Vxj,  Vxfc) 

=  e(Vxj,  Vxfe)  +  (rsrvxj  •  vXfc), 

where  we  have  used  the  periodicity  of  Xk  and  H  in  the  second  equality  and  the  final  equality  follows  from  the 
property  TVxj  =  Vxj.  We  stress  that  T  is  a  projection  on  %  and  therefore  acts  as  the  identity  operation 
on  Li.  Consequently,  for  A  and  S  considered  as  operators  on  Li,  we  have  the  following  weak  equalities 

(c.22)  a  =  rsr  =  sr  =  rs  =  s. 

Since  Vxfc  is  real-valued,  we  have  that  (Vxfc,  Vxj)  =  (Vxj,  Vxfc),  implying  that  S*  is  indeed  a  sym¬ 
metric  matrix.  As  in  Section  C.l,  time  differentiation  and  spatial  convolution  commute,  (— A)_1Ti/>  = 
T(— A)-1,i/>,  for  if  £  J7  [36,  90].  This  and  the  skew-symmetry  of  the  matrix  H  implies  that  S  =  (— A)_1T+H 
is  a  skew-symmetric  operator  on  J- .  Since,  T  is  self-adjoint  on  J-,  TST  is  also  skew-symmetric  on  J- .  Just 
as  in  Section  C.l,  this  implies  that  A*  is  indeed  an  antisymmetric  matrix.  Analogous  to  equation  (2.12)  we 
have  the  following  resolvent  formula  for  Vyj 

(C.23)  Vxj  =  (el  +  A)_1gv,  gJ  =  -THej. 

Equation  (C.23)  follows  from  applying  the  integro-differential  operator  V(A_1)  to  the  cell  problem  in  equa¬ 
tion  (C.18)  written  as  V  •  crVxj  =  — V  ■  H ej,  yielding 

(C.24)  r(el  +  S)VXj  =  -rHep 

The  equivalence  of  equations  (C.23)  and  (C.24)  then  follows  from  the  formula  T V x?  =  Vxj ■  Inserting  the 
resolvent  formula  in  (C.23)  for  Vy,  into  equation  (C.20)  yields  the  following  analogue  of  (2.13) 

(C.25)  S*fe  =  e  {5jk  +  ((el  +  A)-1^.,  (el  +  A)~1gk)) , 

Ajfc  =  (A(el  +  A)  1gJ-,(el  +  A)  1gk), 

We  therefore  have  the  following  corollary  of  Theorem  C.l. 

Corollary  C.3.  The  operator  A  =  TST  with  S  =  (— A)_1T  +  H  displayed  in  equation  (C.20)  is  a 
maximal  (skew-symmetric)  normal  operator  on  the  function  space  T  defined  in  (C.10),  hence  M  =  —  iA  is 
self-adjoint  on  J- .  Let  Q(A)  he  the  resolution  of  the  identity  in  one-to-one  correspondence  with  M.  Define 
the  complex  valued  function  HjkW  =  (Q W9j>9k)>  h  k  =  1,  •  ■  • ,  d,  where  g j  =  —THej  is  defined  in  (C.23) 
and  (•,•)  is  the  LI -inner-product.  Consider  the  positive  measure  gkk  and  the  signed  measures  Regjk  and 
Impjk  associated  with  /ijk(A),  introduced  in  equation  (A. 3).  Then,  for  Uj’Vxj  G  F  and  all  0  <  e  <  oo, 
the  functional  formulas  for  S*k  and  A*k  displayed  in  (C.25)  have  the  Radon-Stieltjes  integral  representations 
displayed  in  equation  (2.14). 

Proof  of  Corollary  C.3.  We  first  establish  that  M  =  —  iA,  with  A  =  TST,  is  a  self-adjoint  operator 
on  T .  Since  T  :  Tip  is  a  projection,  it  acts  as  the  identity  onHx.  We  can  therefore  focus  our  analysis 

on  the  operator  S  =  (— A)-1T  +  H.  As  discussed  above,  since  u  G  T  and  u  =  V  •  H,  the  skew-symmetric 
matrix  H  is  bounded  in  operator  norm  ||H||  <  oo  with  domain  TLp®  Hy,  and  is  therefore  a  maximal  (skew- 
symmetric)  normal  operator,  hence  — *H  is  self-adjoint  on  the  same  domain  [92].  It  is  clear  from  Section  B 
that  — zT  is  self-adjoint  with  domain  Ap.  The  integral  operator  (—A)-1  is  self-adjoint  and  compact  on  "Hy- 
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Since  they  commute  on  Aj-  ®  Hv  [36],  it  follows  that  the  operator  — i(— A)_1T  is  self-adjoint  with  domain 
At®Hv,  hence  (— A)_1T  is  a  maximal  (skew-symmetric)  normal  operator  on  the  same  domain.  Since  — *H 
is  a  self-adjoint  operator  with  domain  Ht  <S>  Hv  and  (—  A)_1T  is  self-adjoint  with  domain  At  <g>  Hv,  the 
operator  S  =  (— A)-1T  +  H  is  self-adjoint  with  domain  D( S)  =  Hr  <8  Hv  D  At  <8  Hv  =  At  <8  Hv  [92]. 
Consequently,  the  operator  M  =  —  lA,  with  A  =  TSr  is  self-adjoint  with  domain  F  =  At  <8  Hx  ■ 

In  Theorem  C.l  we  established  that  the  functions  F( A)  =  (e  +  iA)-1  and  G{ A)  =  *A(e  +  zA)_1  involved  in 
the  functional  formulas  for  S*k  and  A*k  in  equation  (C.25)  are  bounded  for  all  0  <  e  <  oo  so  that  ip  £  $>(F) 
and  <p  £  @(G)  for  all  <p  £  D(M)  when  0  <  e  <  oo.  Since  gj  =  —  THe^  and  1 1 gr ^ 1 1  <  ||r|||[H||  <  oo,  it 
is  clear  that  £  F  C  D(M).  In  the  setting  of  a  mean-zero  velocity  field,  H (t,x)  is  mean-zero  in  time. 
Consequently,  (g.j )  =  0  follows  from  the  Fubini-Tonclli  theorem  [36].  The  conditions  of  the  spectral  theorem 
are  thus  satisfied  which,  from  the  proof  of  Theorem  C.l,  establishes  the  integral  representations  for  S*k  and 
A*fc  displayed  in  equation  (2.14).  From  equation  (A. 6),  the  mass  g{i]k  of  the  measure  g,jk  is  given  by 

(C.26)  9j,gk )  =  (rHej,rHefc>  =  (HrrHei,efc). 

Moreover,  \g°k\  <  ||  H || 2  <  oo  for  all  j,  k  =  1, . . . ,  d.  This  completes  the  proof  of  Corollary  C.3  □. 

We  conclude  this  section  with  the  following  lemma,  which  provides  the  relationship  between  the  effective 
parameter  cr*  defined  in  equation  (C.19)  and  D*  defined  in  (2.9). 

Lemma  C.4.  Let  the  components  D*k  and  a*k,  j,k  =  1, ...  ,d,  of  the  effective  tensors  D*  and  cr*  be 
defined  as  in  equations  (2.3)-(2.9)  and  (C.18)-(C.19),  respectively.  Then  these  effective  tensors  related  by 

(C.27)  cr*  =  [D*]t  +  (H). 

Proof  of  Lemma  C.4.  Recall  the  definition  of  the  function  space  T  =  At  <8  Hx  in  (C.12),  which  we 
will  denote  by  Fx  here.  Denote  by  F,  =  At  <8  H.,  where  H,  C  Hv  is  the  space  of  divergence- free  fields. 
Recall  that  cr  =  el  +  S  where  S  =  (— A)-1T  +  H.  From  equation  (C.18)  we  have  that  Jk  =  crEk  and 
Fk  =  VXfc  +  ek  satisfy  Jk  £  F,  and  Ek  £  Fx.  By  the  mutual  orthogonality  of  the  Hilbert  spaces  Hx  and 
H,  we  have  that  {J j  •  Vy*,)  =  0  for  all  j,k  =  1, . . . ,  d,  which  is  equivalent  to  equation  (C.21).  Consequently, 
from  equation  (C.19)  we  have  (Jj  •  Ek)  =  (Jj  •  ek)  =  cr*k. 

From  S  =  (— A)_1T  +  H  we  have  Se^  =  He^.  By  the  definition  u  =  V  •  H  in  (C.13)  and  the  periodicity 
of  H  and  Xk,  we  also  have  (He^  •  Vyi-)  =  —(ujXk)  via  integration  by  parts.  Therefore,  by  the  skew-symmetry 
of  S  on  Fx,  (VXj)  =  0,  and  equation  (2.9),  we  have 

(C.28)  a*k  =  (J,  •  ek) 

=  ((el  +  S)VXj  •  ek)  +  ((el  +  S)e,  •  ek) 

=  -(Vxi  •  Hefe)  +  ((el  +  H )ej  •  ek) 

=  (wfcXj)  +  e6jk  +  (Hjfe) 

=  D*kj  +  (Hjk), 

which  is  equivalent  to  (C.27).  This  concludes  our  proof  of  Lemma  C.4  □. 

Appendix  D.  An  isometric  correspondence.  A  natural  question  to  ask  is  the  following.  Is  the 
formulation  of  the  effective  parameter  problem  described  in  Theorem  C.l  equivalent  to  that  described  in 
Corollary  C.3?  The  answer  is  in  the  affirmative.  The  correspondence  between  the  two  formulations  is  one 
of  isometry,  and  is  summarized  by  the  following  theorem. 

Theorem  D.l.  The  function  spaces  LF  and  F  defined  in  equations  (C.3)  and  (C.12)  are  in  one-to- 
one  isometric  correspondence.  This  induces  a  one-to-one  isometric  correspondence  between  the  domains 
D(A)  and  D( A)  of  the  operators  A  and  A  defined  in  equations  (2.11)  and  (C.20),  respectively.  Specif¬ 
ically,  for  every  f  £  D(A)  fl  &  we  have  V/  £  D( A)  fl  F  and  ||A/||i  =  ||AV/||,  and  conversely,  for 
each  xp  £  D( A)  n  F  there  exists  unique  f  £  D{A)  D  &  such  that  ip  =  V/  and  ||A ip  ||  =  ||A/||i.  The 
Radon-Stieltjes  measures  underlying  the  integral  representations  of  Theorem  C.l  and  Corollary  C.3  are 
equal,  d(Q(X)gj,  gk)i  =  d(Q(A  )gj,gk),  j,k  =  1  ,...,d,  up  to  null  sets  of  measure  zero,  where  g j  =  Vg^. 
Moreover,  the  operators  A  and  A  are  related  by  AV  =  VA,  which  implies  and  is  implied  by  the  weak 
equality  Q(A)V  =  VQ(A). 
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Proof  of  Theorem  D.l.  We  use  the  formula  u  =  V-  H  displayed  in  equation  (C.13)  and  the  identities 
[V-H  ] •  V/  =  V-[HV/]  and  [V-  =  V-  [He,,]  to  write  the  operator  A  =  (—A)~1(dt  —  u-V)  and  function 

g.j  =  (— A)~1Uj  defined  in  equations  (2.11)  and  (2.12)  as  A  =  (—A)~1(dt  —  V-HV)  and  gj  =  (— A)-1  V- He,,, 
respectively.  Using  the  definition  T  =  —  V(— A)_1V-  and  the  formulas  V(— A)-1c^  =  (— A)_1TV  ,  A  = 
(— A)_1T  +  TH,  and  g  •  =  — T He,,,  displayed  in  equations,  (C.17),  (C.22),  and  (C.23),  respectively,  we  have 

(D.l)  VA=[(- A)-1T  +  TH]V  =  AV,  V9j=gy 

Consequently,  by  applying  the  differential  operator  V  to  both  sides  of  the  formula  (e  +  A)xj  =  gj  of  (2.12), 
we  obtain  the  formula  (el  +  A )V\j  =  <7j  of  equation  (C.23). 

Since  the  function  spaces  &  and  T  differ  only  in  the  characterization  of  the  spatial  variable  x,  we  now 
discuss  the  relationship  between  the  Hilbert  spaces  My  and  defined  in  equations  (C.2)  and  (C.ll), 
respectively.  We  will  temporarily  denote  the  inner-product  induced  norms  of  these  Hilbert  spaces  by  ||  •  ||i 
and  ||  •  ||,  respectively.  The  Sobolev  space  My  in  (C.2)  is  the  closure  in  the  norm  ||  •  ||i  of  the  space 
C2(V)  of  all  twice  continuously  differentiable  periodic  functions  in  My  [14].  Functions  in  My  need  not  be 
differentiable  in  the  classical  sense.  Instead,  /  G  My  has  derivatives  dkf/dxk  G  L2(V)  for  k  =  1,2,  defined 
by  dkf/dxk  =  Hindoo  dkfn/dxk  where  /„  G  C2(V)  are  Cauchy  in  the  norm  ||  •  ||i  and  converge  to  /  in 
i2(V)  [63], 

Let  {/„}  be  such  a  Cauchy  sequence,  i.e. ,  for  all  e  >  0,  there  exists  an  IV  G  N  such  that  for  all  n,m  >  N 
we  have  that  ||/ra  — /m||i  <  e  and  lim^oo  /„  =  /,  where  /  G  My.  Since  /„  G  C2(V)  and  /„  is  also  V-periodic, 
we  have  that  A-1  A /„  =  /„  [63].  Consequently, 

(D.2)  e  >  ||/„  -  fm\\i  =  ||  VA-1A/„  -  V/m||  =  ||TV/„  -  V/m||, 

which  implies  that  TV/  =  V/.  Moreover,  ||V/||2  =  (V/  •  V/)  =  ||/||2  <  oo.  Consequently,  for  every 
/  G  My  we  have  V/  G  'Hx-  Conversely,  ip  G  Hx  implies  that  ip  =  Tip  =  V/,  where  we  have  defined  the 
scalar- valued  function  /  =  A  ~XV -ip .  Since  if)  =  V/,  the  My  norm  of  /  satisfies  ||/||2  =  ('»/>•'!/’)  =  ||'0||2  <  oo 
so  that  /  G  My.  Moreover,  /  is  uniquely  determined  by  il>  (up  to  a  zero  Lebesgue  measure  equivalence 
class),  since  if  f\  =  A_1V  •  ip  and  fi  =  A_1V  •  ip  then  Tip  =  ip  implies  that  ||/i  —  /2II1  =  \\ip  —  ip\\  =  0. 
Consequently,  for  every  ip  G  Hx  there  exists  unique  /  G  My  such  that  ip  =  V/.  In  summary,  the  Hilbert 
spaces  My  and  are  in  one-to-one  isometric  correspondence,  which  we  denote  by  My  ~  Hx-  This,  in 
turn,  implies  &  ~  J- . 

We  now  return  to  our  previous  notation,  where  ||  •  ||i  and  ||  •  ||  denote  the  norm  induced  by  the  M- 
inner-product  and  H-inner-product,  respectively.  We  demonstrate  that  the  one-to-one  isometry  between 
&  and  T  induces  a  one-to-one  isometry  between  the  domains  D(A)  and  D( A)  of  the  operators  A  and  A. 
This,  in  turn,  follows  from  the  one-to-one  correspondence  between  a  self-adjoint  operator  and  its  resolution 
of  the  identity  discussed  in  Section  A,  leading  to  equation  (A. 9).  More  specifically,  the  domain  D{M)  of  the 
self-adjoint  operator  M,  for  example,  comprises  those  and  only  those  elements  /  of  M  such  that  the  Stieltjes 
integral  J  A2  d||Q(A)/||2  is  convergent,  and  when  /  G  D(M)  the  element  Mf  is  determined  by  the  relations 
in  equation  (A. 9).  Since  A  =  iM  it  is  clear  that  D(A)  =  D(M).  We  already  established  in  Section  C.l  that 
&  C  D(A),  and  similarly  for  the  operator  A. 

Let  /  G  D(A)  (~l  & .  From  the  relation  &  ~  J7,  we  have  that  V/  G  T,  so  from  equation  (D.l) 

(D.3)  \\Af\W  =  (Af,  Af h  =  (VAf  ■  VAf)  =  (AV/  •  AV/)  =  || AV/||2. 

Consequently,  from  equation  (A. 9)  we  have  that 

(D.4)  J  A2d||Q(A)/||2  =  J  A2  d||Q(A)V /||2, 

and  the  convergence  of  the  integral  on  the  left-hand-side  of  (D.4)  implies  the  convergence  of  the  integral  on 
the  right-hand-side  which,  in  turn,  implies  that  V/  G  D( A). 

Conversely,  let  ip  G  D( A)  D T .  From  the  relation  &  ~  J7,  there  exists  unique  /  G  &  such  that  i/j  =  V/, 
and  equation  (D.l)  then  implies  that 


(D.5) 


II  AV>  ||2  =  (AV/,  AV/)  =  (VAf,  VAf)  =  (Af,  Af),  =  \\Af\\2. 
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Again,  equation  (A. 9)  implies  that  (D.4)  holds,  and  the  convergence  of  the  integral  on  the  right-hand-side 
of  (D.4)  implies  the  convergence  of  the  integral  on  the  left-hand-side  which,  in  turn,  implies  that  /  £  D(A). 

In  summary,  for  every  /  £  D{A)C\£F  we  have  V/  £  D( A)  and  ||A/||2  =  ||AV/||2.  Conversely,  for  every 
xf  £  D( A)  fl  T,  there  exists  unique  /  £  D(A)  such  that  xf  =  V/  and  ||  Axp  ||2  =  ||A/||2.  This  generates  a 
one-to-one  isometric  correspondence  between  the  domains  D( A)  and  D{A). 

We  now  show  that  this  result  implies,  and  is  implied  by  the  weak  equality  VQ( A)  =  Q(A)V,  where 
Q{ A)  and  Q(A)  are  the  resolutions  of  the  identity  in  one-to-one  correspondence  with  the  operators  A  and 
A,  respectively.  From  equation  (D.4)  and  the  linearity  properties  of  Radon-Stieltjes  integrals  [92],  we  have 
that 

/oo  poo 

A2  d(||Q(A)/||?  -  ||Q(A)V/||2)  -  /  A2  d«[VQ(A)  -  Q(A)V]/  •  V/)). 

-oo  J  —  oo 

Equation  (D.6)  implies  that  for  all  /  £  D(A)f\&  <*=>•  V/  £  D(A)flJr  we  have  d||Q(A)/||2  =  d||Q(A)V/||2, 
up  to  sets  of  measure  zero.  Moreover,  the  equality  VQ(A)  =  Q(A)  V  holds  in  this  weak  sense.  Conversely, 
assume  that  Q( A)  and  Q(A)  are  the  resolutions  of  the  identity  in  one-to-one  correspondence  with  the  operators 
A  and  A  and  that  VQ(A)/  =  Q(A)V/  for  every  /  £  D{A)C\&  -<=>•  V/  £  D(A)nJr.  Then  equation  (D.6) 
holds  and  implies  equation  (D.4).  Equation  (A. 9)  then  implies  that  ||AV/||2  =  ||A/||2  =  ||VA/||2,  which 
implies  that  AV  =  VA  in  this  weak  sense.  Since  gk  £  D(A)  and  gk  £  D( A)  with  gk  =  "Vgk ,  this  result 
implies  that  the  Radon-Stieltjes  measures  underlying  the  integral  representations  of  Theorem  C.l  are  equal 
to  that  of  Corollary  C.3  d||<3(A)<?fc||i  =  d||Q(A)gfc||  up  to  null  sets  of  measure  zero,  for  all  j,k  =  1 ,d. 
This  concludes  our  proof  of  Theorem  D.l  □. 

D.l.  Discrete  integral  representations  by  eigenfunction  expansion.  The  integral  representa¬ 
tions  of  Theorem  C.l  and  Corollary  C.3  displayed  in  equation  (2.14),  involve  a  spectral  measure  (ijk-, 
j,k  =  1 , . . . ,  d,  which  has  discrete  and  continuous  components  [84,  92].  In  this  section,  we  review  these 
properties  of  p,jk  and  provide  an  explicit  formula  for  its  discrete  component.  Towards  this  goal,  we  summa¬ 
rize  some  general  spectral  properties  of  the  self-adjoint  operators  M  =  —iA  and  M  =  —  iA  on  the  function 
spaces  &  and  J r,  which  are  densely  defined  on  the  associated  Hilbert  spaces  Aff  and  TL,  given  in  equa¬ 
tions  (C.3)  and  (C.12),  respectively.  We  will  focus  on  the  operator  M  and  the  Hilbert  space  Atf,  as  the 
discussion  regarding  M  and  TL  is  analogous. 

Recall  from  equation  (A. 9)  that  the  domain  D(M )  of  the  self-adjoint  operator  M  comprises  those  and 
only  those  elements  /  £  Aif  such  that  ||M/||2  =  A2  d||Q(A)/||f  <  oo,  where  Q( A)  is  the  resolution  of  the 
identity  in  one-to-one  correspondence  with  M  [92],  The  integration  is  over  the  spectrum  E  of  M,  which  has 
continuous  Econt  and  discrete  (pure-point)  Epp  components,  E  =  Econt  U  Epp  [84,  92].  We  first  focus  on  the 
discrete  spectrum  Epp. 

The  /  £  ^  0,  satisfying  Mf  =  A /  with  A  £  Epp  are  called  eigenfunctions  and  A  is  the  corresponding 

eigenvalue.  Since  M  is  self-adjoint,  A  is  real- valued  [92].  The  span  of  all  eigenfunctions  is  a  countable 
subspace  of  Jif  [92].  Accordingly,  we  will  denote  the  eigenfunctions  by  ipi,  l  =  1, 2,  3, . . .,  with  corresponding 
eigenvalues  A;.  Eigenfunctions  corresponding  to  distinct  eigenvalues  are  orthogonal  and  can  be  normalized 
to  be  orthonormal  [92],  i.e.  if  Mpi  =  A  ipi  and  Mpm  =  \mpm  for  A  i  Am,  then  (pm,pn)  l  =  5mn.  There 
can  be  more  than  one  eigenfunction  associated  with  a  particular  eigenvalue.  However,  they  are  linearly 
independent  and,  without  loss  of  generality,  can  be  taken  to  be  orthonormal  [92].  Consequently,  associated 
with  each  eigenfunction  pi  is  a  closed  linear  manifold,  which  we  denote  by  A 4{pi)-  When  l  p  m,  A i(pi)  and 
M(tpm)  are  mutually  orthogonal.  Set  £  =  ®f£iM(pi)7  M  =  £  ©  {0},  and  let  J\f  =  M ^  be  the  orthogonal 
complement  of  A4  in  Af.  All  the  properties  of  A4  and  M  that  are  relevant  here  have  been  collected  in 
the  following  theorem  [92],  which  provides  a  natural  decomposition  of  the  Hilbert  space  Atf  in  terms  of  the 
mutually  orthogonal,  closed  linear  manifolds  M  and  A A,  and  leads  to  a  decomposition  of  the  measure  gkk 
into  its  discrete  and  continuous  components. 

Theorem  D.2  ([92]  pages  189  and  247).  One  of  the  three  cases  must  occur: 

1.  £  =  0  and  A4  =  {0}  has  dimension  zero;  M  =  At?  has  countably  infinite  dimension.  There  exists 
an  orthonormal  set  {xpm},  m  =  1,2,3,...,  and  mutually  orthogonal,  closed  linear  manifolds 
which  determine  A f  according  to  Af  =  ®^=1A 

2.  £  contains  an  incomplete  orthonormal  set  {pi}  so  that  both  A4  and  Af  are  proper  subsets  of  Aif , 
Af  having  countably  infinite  dimension  and  A4  having  finite  or  countably  infinite  dimension.  There 
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exists  an  orthonormal  set  {ipm}  in  Af .  The  closed  linear  manifolds  A ’i(ipi)  and  Af  {ipm)  are  mutually 
orthogonal  and  together  determine  Iff  according  to 

W  =  ©~=1W(V>m),  = 

3.  £  contains  a  complete  orthonormal  set  {tpi};  AA  =  Iff  has  countably  infinite  dimension;  A f  =  {0} 
has  zero  dimension.  In  this  case ,  the  closed  linear  manifolds  AA(ipi)  are  mutually  orthogonal  and 
together  determine  AA  according  to  AA  =  ®}fi1AA(ipi). 

In  each  of  these  three  cases ,  the  closed  linear  manifolds  AA  and  Af  reduce  M,  i.e.,  M  leaves  both  AA  and 
Af  invariant  in  the  sense  that  if  f  £  D(M )  and  f  £  Af  then  M f  £  Af,  and  similarly  for  AA.  In  cases  (2) 
and  (3),  a  necessary  and  sufficient  condition  that  an  element  ipi  £  Iff  be  an  eigenfunction  with  eigenvalue 
X i,  is  that  the  function  ||Q(A)</3/||j  is  constant  on  each  of  the  intervals  —  oo  <  A  <  A i  and  Xi  <  A  <  oo  [92]. 
Moreover,  a  necessary  and  sufficient  condition  that  f  £  AA,  /  ^  0,  is 

OO  OO 

(D.7)  /  =  £</,W>  i<Pi,  ll/ll?  =  El</’^)i|V°, 

;= i  ;= i 

and  similarly  for  f  £  Af  with  orthonormal  set  {ipm}-  In  cases  (1)  and  (2),  a  necessary  and  sufficient 
condition  that  ip  0  be  an  element  of  Af  is  that  ||Q(A)^||j  be  a  continuous  function  of  X  not  identically 
zero  [92]. 

Let  f  be  an  arbitrary  element  of  Iff ,  and  g  and  h  be  its  (unique)  projections  on  A4  and  Af ,  respectively, 
then  the  equation 

(D-8)  ll<9(A)/||?  =  HQ(A)gili  +  ||Q(A)/i||i, 

d||Q(A)/||?  =  d||Q(A)ff||?  +  d||Q(A)/l||? 

is  valid  and  provides  the  standard  resolution  of  the  monotone  function  ||Q(A)/||f  into  its  discontinuous  and 
continuous  monotone  components,  as  well  as  the  decomposition  of  the  measure  d||(5(A)/||^  into  its  discrete 
and  continuous  components. 

We  now  use  the  mathematical  framework  summarized  in  Theorem  D.2  to  provide  explicit  formulas  for 
the  discrete  parts  of  the  integral  representations  for  5*k  and  A*fc,  displayed  in  equation  (2.14).  Recall  the 
cell  problem  in  equation  (2.8)  written  as  in  (C.4),  (e  +  A)xj  =  gj ■  Here  A  =  iM  is  defined  in  (2.11), 
gj  =  (— A)"1^,  and  Uj  is  the  jth  component  of  the  velocity  field  u,  j  =  1  ,...,d.  Moreover,  we  have 
Xj,gj  £  C  Iff  and  &  C  D(A).  We  stress  that  the  arguments  presented  here  are  more  subtle  than  those 
typically  used  for  bounded  operators  in  Hilbert  space.  The  reason  being  is  that  a  bounded  linear  operator 
commutes  with  all  the  infinite  sums  encountered  here,  by  the  dominated  convergence  theorem  [36].  However, 
for  the  operator  A ,  we  must  instead  rely  on  general  principles  of  unbounded  linear  operators  in  Hilbert  space. 

Let  Xj  and  xf  be  the  (unique)  projections  of  Xj  011  -Ad  and  Af,  respectively,  with  Xj  =  Xj  +  xf  and 
similarly  for  gj.  Since  A  =  iM  is  a  linear  operator,  we  have  that  Axj  =  Axj  +  A\f  ■  From  Theorem  D.2,  the 
linear  manifolds  A4  and  Af  both  reduce  A,  which  implies  that  Axj  £  AA  and  A\f  £  Af.  From  equation  (D.7) 
we  then  have  A\j  =  Y^i(AXji  Tl)  i  fi  and 

(D-9)  Xj  =  ^2(Xj,Ti) i  Ti  +  xf ,  AXj  =  lMXj,Vl) l  *Pi  +  Axf 

i  i 

where  we  have  used  that  (Axj,ipi) i  =  —(xj,Aipi) ±  =  — ( Xj,iXnpi)\  =  iXi(xj,ipi)\.  From  the  cell  problem 
(e  +  A)xj  =  gj  we  therefore  have 

(D.10)  e^ixj^tim  +^2iXi(xj,<pi)m  +  {e  +  A)xf  =9j  +gf, 

i  i 

where  (e  +  A)xf,gf  £  Af.  Of  course,  each  f  £  Af  can  be  represented  [92]  as  /  =  Ylmif^^ifpm,  where 
{ipm}  is  the  orthonormal  set  defined  in  Theorem  D.2,  though  we  have  suppressed  this  notation  in  the  above 
equations  for  simplicity.  By  the  mutual  orthogonality  of  the  linear  manifolds  AA  and  Af,  the  completeness 
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of  the  set  {</?;}  U  {ipm}i  and  the  Parseval  identity,  taking  the  inner-product  of  both  sides  of  equation  (D.10) 
with  ipn  yields 


(D-ll) 


{Xj,fn)  1 


(dj  I  yra)l 


0  <  £  <  OO. 


Recall  the  representations  S*k  =  e(Sjk  +  {XjiXk) l)  and  A*k  =  (A\j,Xk) l,  j,k  =  1  displayed  in 

equation  (2.11).  Writing  \j  =  Xj  +  xf  and  A\j  =  A\j  +  A\f,  the  mutual  orthogonality  of  the  linear 
manifolds  M  and  Af,  which  both  reduce  A ,  implies  that  {XjiXk) l  =  (XjiXk) l  +  (xft  Xk  )i  and  {A\j,Xk)i  = 
(Axj,Xk)i  +  (Axf,Xfc)i-  Consequently,  from  equations  (D.9)  and  (D.ll),  the  completeness  of  the  set 
{c,Ci}  U  {ipm},  and  the  Parseval  identity,  we  have 


(D.12) 


(Xj,Xk) l  -  ixf,xi)  1  =  ^2(Xj,<Pl)i(Xk,<Pl)  1  =  X! 


(gj,<pi)i(gk,vi)1 

e2  +  A? 


-  (Axf  ,Xfe)i  =  53  A 


A  {gj,g’i)i{gk,‘Pi}1 


£2  +  \? 

i  i  1 

Since  Xj  and  A\j  are  real- valued,  just  as  in  equation  (A. 8),  we  have  that 


(D.13) 


(Xj,Xfc)i 


(Axj,Xk) i  - 


\Aj?Afc/l  £.2  _j_  ^2 


X^d-Xj  ,Xfc)l  r2  _|_ 


The  right  hand  sides  of  the  formulas  in  equation  (D.13)  are  Radon-Stieltjes  integrals  associated  with  a 
discrete  measure.  The  terms  {xfiXk)i  and  {A\f  ,Xk) i  a^so  have  Radon-Stieltjes  integral  representations 
involving  the  continuous  measure  d(Q(X)gj- ,  gk)i  via  equation  (2.14).  We  note  that  from  the  decomposition 
9j  =  9j  +  9 we  have  (§j,<pi) i  =  {gj,fi) i-  A  useful  property  of  the  inner-product  ( gj,<Pi)i  and  the  form  of 
gj  =  (—A)~1Uj  is  that  (see  equation  (C.5)) 

(D.14)  (gj,Vi)i  =  {Uj,<pi)2- 

This  property  is  used  in  Section  3  to  calculate  S*fc  and  A*fc  for  a  large  class  of  velocity  fields. 
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